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I. 



INTRODUCTION 



Although combat is a complex random process, analysts 
frequently model it with deterministic Lanches ter- type 
equations for reasons of mathematical tractability and com- 
putational convenience. Such deterministic models are being 
used extensively today in defense planning studies. However, 
one can include the various random aspects of combat (e.g., 
uncertainty in enemy force level, random occurrence of 
casualties, etc.) in models and develop so-called stochastic 
models of the combat process. Thus, a basic way of classi- 
fying a combat model is whether the model is a deterministic 
model or a stochastic one. 

Moreover, there are fundamental differences between 
deterministic and stochastic models. Deterministic dynamic 
combat models predict the future with certainty; for given 
initial conditions there is no question about what the state 
of the conflict will be at any future time. On the other 
hand, stochastic dynamic models only tell one the "chances" 
of what will happen. With a stochastic model one does not 
know the future with certainty. That is, it cannot be 
guaranteed that a specific state will be reached in the 
future — there is only a "probability" it will be reached. 

Each type of model has inherent advantages and disadvantages 
which are partially summarized in Table I. 



12 



Deterministic Models 



Advantages Disadvantages 

1. Requires one solution of the 1. A further abstraction 

model for a useful outcome from reality 

2. Less expensive to run 

3. Interactions are more 
easily analyzed 



Stochastic Models 



Advantages 



Disadvantages 



1. Closer to reality 



1. Requires many replica- 
tions to provide 
statistically 
significant results 

2. More expensive to run 

3. Interactions are much 
more difficult to 
analyze 



TABLE I. Advantages and Disadvantages 

of Deterministic and Stochastic 
Models 
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The solutions to stochastic models quickly become quite com- 
plex as the number of combatants increases, and are not very 
useful; so one must use other methods to obtain insights into 
the dynamics of combat. The two most commonly used methods, 
Monte Carlo simulation and finite difference approximation, 
both possess certain disadvantages. Monte Carlo combat simu- 
lations require that the battle be replicated many times in 
order to obtain a good statistical estimate of the probable 

cc 

course of battle, while finite different approximation methods 
usually require digital computer implementation (including 
development of a computer program) . 

On the other hand, solutions to deterministic models can 
sometimes be useful for developing insights into the dynamics 
of combat. Even when they are not, the model must be solved 
only once by finite difference approximation as opposed to 
the many replications required for Monte Carlo simulation. 

Also, finite difference approximations of deterministic models 
are generally less costly in terms of time, and manhours 
required to implement the computer program, than the finite 
difference approximations of stochastic models. The economic 
costs of using the models then, is the basic reason most models 
of combat used in the past have been deterministic. But can 
the complex random process of combat be adequately represented 
by a deterministic model? Does one destroy any of the essential 
features of the random combat process by considering a deter- 
ministic model as representing the "average" course of combat? 
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The objective of this thesis is to consider this question by 
investigating the differences in results obtained from a 
deterministic and a stochastic model of combat. 

The importance of this objective stems from the fact that 
both types of models are used extensively today in the decision 
making process in the U.S. Army. For example, the Dynamic 
Tactical Simulation, DYNTACS-X, and the Battalion Differential 
Model, BLDM, widely known as the Bonder/IUA model, are two 
models that simulate battalion level mid-intensity armored 
combat. DYNTACS is an event sequenced stochastic model that 
is being or has been used by the Rock Island Arsenal to evalu- 
ate proposed mobility improvements to the M60 tank, the family 
of scatterable mines concept, cannon launched guided projec- 
tiles, remotely piloted vehicles, and the XM1 tank. Bonder/ 

IUA is a deterministic Lanchester-type model that is being or 
has been used by the Rock Island Arsenal in sensitivity analyses 
and evaluation of the MBT70 study, several anti-armor automatic 
cannon concepts, and the Low Dispersion Automatic Cannon study. 

A study done by the Rock Island Arsenal [Ref. 5] shows that 
there are differences in the results of the two models with 
"equivalent" inputs. With the present state of the art, it 
is not possible to predict with any accuracy when significant 
differences between the results of the models will occur. It 
is felt that this thesis may provide trends which will allow 
a user to gain insights into these differences and more 
accurately determine when a deterministic model represents 
combat "satisfactorily". 
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Models have had wide application in the military decision 
making process, and the realm of possible future applications 
is even larger. For example, the Engineer Strategic Studies 
Group [Ref. 10] has, as its name implies, evaluated several 
aspects of military engineering, but with little emphasis on 
the use of combat models . Much of this has been caused by a 
dearth of models that could be applied to engineer studies. 
However, this is being altered; for example, the U.S. Army 
Concepts Analysis Agency [Ref. 30] has accepted as a future 
tool a model developed by ESSG, the Fortification and Obstacle 
Effect Simulation, 1975. This is a high resolution model 
used in analyzing obstacle effectiveness. The model is part 
of a research effort by ESSG to design a method for assessing 
the contribution that obstacles make to combat operations. 

This model was developed using ideas from other models, but 
it was based on the identification of specific objectives. 

From these objectives measures of effectiveness were developed, 
and from these measures the model was developed. It is possi- 
ble, with the identification of proper objectives, to develop 
models for other phases of military engineering and, more 
generally, for the evaluation of other phases of combat. It 
is therefore felt that combat models do, and will continue to 
have an ever increasing application to the evaluation of 
military systems. 

As was noted, however, it is not always clear which type 
of model should be applied — a stochastic or a deterministic 
model. From an economic standpoint, the advantages of the 
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deterministic model weight heavily in its favor. If the 
deterministic results do not significantly differ from those 
of the stochastic model, it seems the former should be used. 
The question is, when are the models "similar"? 
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II. LITERATURE SURVEY 



This chapter briefly reviews Lanchester's well-known 
equations of combat and some simple extensions, with the 
discussion focused on stochastic formulations of the combat 
attrition process. The review consists of three parts; 

(A) Deterministic Models, (B) Stochastic Models, and 
(C) Comparisons of the two types. 

A. DETERMINISTIC MODELS 

In 1914 Lanchester [Ref. 19] hypothesized that under 
conditions of "modern warfare" combat between two homogeneous 
forces could be described by the equations 



dx 

dt 



ay 



( 1 ) 



dy = 

dt 



bx 



( 2 ) 



with initial conditions 



x(t = 0) = x Q (3) 

y(t = 0) = y Q (4) 



where a and b are commonly referred to as the Lanchester 
attrition rate coefficients and x(t) and y(t) are force 
levels. One set of assumptions that has been hypothesized 
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(proposed by H. Weiss [Ref. 35] and cited by Dolansky 

"7 

[Ref. 9]) to yield equations (1) and (2) is 

1. Two forces are engaged in a battle. The units 
are homogeneous , but the rate of attrition may 
be different for each force. 

2. Each firer on either side is within the range 
of all targets on the other side. 

3. Attrition rate coefficients are constant. 

4 (a) . Firers have perfect knowledge of the target 
locations and disposition so that they fire only 
at live targets, and fire can be shifted 
instantaneously. 

5(a). Fire is uniformly distributed over remaining 
targets . 

The state solution (time "independent" solution relating 
the force levels ) is given by 



b(x Q 2 - x(t) 2 ) = a(y Q 2 - y(t) 2 ) (5) 



hence the name "square-law" attrition process. 

The time solution for equations (1) and (2) is given 
by (see Morse and Kimball [Ref. 22]) 



x(t) = x Q cosh (-y ab t) - ■y' a/b y Q sinh (y/ ab t) (6) 



y (t) = y cosh ( y/ ab t) - y b/a x Q sinh (y/ ab t) (7) 



The square-law was developed as a by-product of the 
original work by Lanchester to quantitatively justify the 
principle of concentration under certain combat conditions. 
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Lanchester was contrasting these conditions with the condi- 
tions of "ancient warfare" between two homogeneous forces. 
He postulated that the latter could be described by the 
equations 



dx 

at ’ ‘ axy 


(8) 


I - - 


(9) 


with initial conditions 




x ( t = 0) = x Q 


(10) 


y ( t = 0) = y Q 


(11) 


One set of assumptions that has 
proposed by H. Weiss [Ref. 3ff] 


been hypothesized (also 
and cited by Dolansky [Ref. 



to yield equations (10) and (11) are very similar to the 
assumptions listed previously for the square-law. The first 
three assumptions are identical to the first three square- 
law assumptions. 

4(b). Each firer knows the area that contains 
targets, but does not know exact target 
locations or the consequences of his fire. 

5(b). Fire from surviving units is distributed 
uniformly over the target area. 

The state solution is given by 



b(x Q - x ( t ) ) = a(y Q - y(t)) (12) 
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hence the name "linear-law" attrition process. 



The time solution for equations (8) and (9) is given 
by (see G. Weiss (Ref. 35]) 

(y “b/ax ) x exp( (b/ax -y ) at) 

~ (y Q -b/ax o exp ( (b/ax Q -y o ) at) 

y(t) = y o ~b/ax o + b/ax(t) (14) 



Others have suggested general forms of homogeneous force 



K 

models given by 
dx 



dt 



F ( t, x,y ) 



= 

dt 



G( t,x,y) 



(15) 



(16) 



For example, Willard [Ref. 37] considered the following 
Lanchester-type equations: 



— = - ax c f ( t,x,y) (17) 

ft = - by C f (t,x,y) (18) 



Many "general" forms contain variable coefficients (e.g., 
a = a(t) and b = b(t)). Taylor [Ref. 30] has shown that 
except for some very special cases, the solutions to varia- 
ble coefficient Lanchester-type equations are very complex 
and in many cases of little practical use. 
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B. STOCHASTIC MODELS 



Koopraan [Ref. 22] suggested a reformulation of the 
attrition process in stochastic form. The resulting sto- 
chastic attrition process has been appropriately called by 
him the Lanchester stochastic process . Others have subse- 
quently employed a stochastic analysis of combat. With the 
inclusion of random variations in the attrition process, 
analysts have attempted to better represent the complex 
random process of combat with the goal of gaining insights 

into combat not available from deterministic models. 

w 

Snow [Ref. 27] showed that the following assumptions 
yield the square-law attrition process: 



(Al) P (t+At,m, n: t ,m+l ,n) = anAt 

(A2) P ( t+At,m,n: t,m,n+l) = bmAt 

(A3) Probability of more than one casualty 
occurring in a time of length At is of 
the order of magnitude o(h) where 

lim ■ ■ = 0. 

h-0 h 



From these, he developed the Chapman-Kolmogorov forward 
equations : 



dP (t,m,n) 
dt 



anP(t,m+l,n) + bmP(t,m,n+l) 



(an+bm) P ( t ,m, n) 
(19) 



If one assumes, instead 

(Al) P (t+At,m,n: t,m+l,n) = F(t,m+l,n)At 
(A2) P ( t+At ,m, n : t ,m, n+1) = G(t,m,n+l)At 
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(A3) Probability of more than one casualty 
occurring in a time of length At is of 
magnitude o(h) . 



then the resulting Chapman-Kolmogorov forward equations are: 



In this thesis, this process will be called the general 
Lanches ter- type stochastic attrition process . 

Work has been done by several authors to obtain general 
solutions for specific stochastic attrition processes. 

Brown [Ref. 7] attempted to find the time state probabili- 
ties (P(t,m,n)) for a class of homogeneous force stationary 
Markov attritions processes (that is, when F(t,m,n) = A(m,n) 
and G(t,m,n) = B(m,n)). He concluded that the expression 
derived is too complicated for practical use. Isbell and 
Marlow [Ref. 14] developed a general solution to (20) for 
a square-law stochastic attrition process. If 
F(t,m,n) = an + 3m and G(t,m,n) = bm + an (Isbell and Marlow 
referred to a and 3 as operational loss rates) , with 
a + a = b + 3 / then 



dP (t ,m, n) 
dt 



= F(t,m+l,n) P (t,m+l,n) + G (t,m,n+l) P (t,m, n+1) 



(F(t,m,n) + G(t,m,n) ) P (t,m,n) 



( 20 ) 




m +n -m-n 
o o 



(21) 



• exp (- (b+3) (m o +n Q ) t) 
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where F (m, n :m Q , n Q ) is developed from a recursive relationship 



m +an 

F(m,n:m ,n ) = yr— . ' ° --- , ° F(m,n:m -l,n ) 

o o (b+p)m +(a+a)n o o 

o o 



n +bm 

+ (bH-e)n.° + (a°a)n 0 F<».n!m 0 ,n 0 -1) (22) 



where F (m o ,n o :m o ,n o ) = 1. 

Clark [Ref. 8] developed a general solution for the time 
state probabilities for the linear-law stochastic attrition 
process (i.e., F(t,m,n) = amn and G(t,m,n) = bmn) ; 



m 



, , . m -m n -n 

n '/ix k-n+i-m o . o 
o o (-1) J a b 



P(t,m,n) = Z Z ( 
j=m k=n 



(m o )! <n o ) 1 



l 



n -n+m -m 

(a+b) ° ° m! n! (k-n) ! (j-m) ! (n Q -k) ! 



m -3 

ii . {j-g- . f . ; n (n 0 -3 + ^) 

£ =1 U-k t) l=1 ( j +k+ ^) 



exp (- (a+b) jkt 
( 23 ) 



(However, he did not compare his results with those for the 
probability of winning (i.e., "true" absorption probabili- 
ties) . This expression is too complicated to provide any 
insights by itself. This author has no knowledge of any 
other general solutions developed for the time state proba- 
bilities. It is, however, a straight forward task to compute 
P(t,m Q ,n) and P(t,m,n Q ), which are used in the included 
computer program as a partial check for accuracy. 
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Work has also been done to develop exact expressions for 
the probability of winning, for a stationary Markov attrition 
process. The exact expression for the probability of winning 
for the linear-law stochastic attrition process is easily 
obtained from classical random walk results [Ref. 28]: 



n n 

o c 

P v = Z P(o,n:m ,n ) = £ 

x n=l ° ° 



m +n -n-1 

( ° ° ) 

, 1 m -1 ' 'a+b 

n=l o 



m n -n 

o.o 

) (— — ) 

1 ' a+b 

(24) 



Smith [Ref. 26] extended work done by Brown [Ref. 7] and 
developed the explicit expression for the square-law stochas- 
tic attrition process: 



m 



P v = Z P(m,0:m ,n ) 
x , o o 

m=l 



(25) 



where 



P(m,0:m o ,n o ) 



bm 

[ — 

an +bm 
o o 



o 

] Z [ 
j=m 



(- 1 ) 



m -i 
o J 



T (m - j + 1) 
o J 



m +n -m- 1 

j ° r (b/a j+1) . 

r(n Q +b/a j+1) (j-m+1) 



(26) 



This is an impressive looking expression, but from a practical 
viewpoint, it is not very useful. A computer is needed to 
solve this equation for any appreciable force sizes, but 
one immediately runs into overflow and underflow problems. 

To solve these complications, recursive relationships must 
be used; but if one looks carefully at the relationships . 
used, he is lead back to the initial recursive relationships 
used to obtain the solution: 
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+ 



( 27 ) 



P (0 ,n:m Q ,n o ) 



A(m ,n ) 
o o 



P (0 / n:m o -l,n o ) 



A(m o' n o )+B(m o' n o ) 



B (m o , n 0 ) 



P (0 ,n:m ,n -1) 
o o 



A (m o ' n o )+B(m o' n o ) 



It was found by this author that the most efficient method 
to solve the stochastic model for the probability of winning 
was by using the above recursive relationship, or one devel- 
oped by Springall [Ref. 28], not the general force level 
solutions (See Appendix A) . 

Springall developed a recursive relationship . for the 
distribution of the time of battle termination: 



r ° 

E(t ) = l 





(28) 



where 




A (m o , n o )+B(m o' n o ) 



1 



(29) 





(30) 



A(m Q ,n) +B (m Q ,n) 





(31) 



m,n 



A (m, n Q ) +B (m,n Q ) 
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X 



A(m+l,n)X^T +B ( m , n+ 1 ) X ^ ^ . 

(0) _ m+1 , n ' m,n+l 

~m, n A(m, n) +B (m, n) 



( 32 ) 



for i > 1 



iX 



(i-1) 



X 



(i) 



m , n 
o' o 



m ,n 

o o 

A (m , n ) +B (m , n ) 
o o ' o o 



( 33 ) 



X 



(i) 



m o' n 



iX (l_1) +B(m ,n+l)X (l) , , 

m ,n o m ,n+l 

o o 

A(m ,n)+B(m ,n) 



( 34 ) 



X 



(i) 



iX (l_1) +A(ra+l,n )X^J 
m,n o m+l,n 

o o 



m, n A(m, n ) +B (m, n ) 



( 35 ) 



( (i) 
m f n 



iX (l 1) +A(m+l,n) X ( i! +B (m, n+1) X (l) 
m f n ' m+ 1 , n m, n+1 

A (m, n) +B (m, n) 



( 36 ) 



From studying the equations presented in this section, one 
gains an appreciation for the increase in complexity of the 
model by the inclusion of uncertainty in the attrition 
process. It is apparent that, although exact solutions may 
be obtained for some stochastic models, numerical results 
are not readily generated by hand. A computer is therefore 
essential for the numerical solution of the simple stochastic 
models . 

It is surprising that more use has not been made of the 
computer in investigating stochastic models, especially the 
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use of computer graphics. For example, the five following 
figures are plots of the time state probabilities for a 
stochastic model at 3, 25, 50, 75, and 100% of the time to 
completion of the equivalent deterministic model. It is 
believed that something of this nature has not been utilized 
in the open literature before. It gives some interesting 
insights into the progress of a stochastic model. The state 
probability seems to start as a spike at (m o ,n Q ) and pro- 
gress^over time somewhat as a drop of water would. As time 
increases, the drop spreads and moves toward the banks 
(boundaries of the state space) . The angle at which the 
wave moves is determined by the relationship between the 
attrition rates of the process (in this case, the model is 
the square-law stochastic attrition process, so the rela- 
tionship is between a, b, m, and n) . As the "wave" hits 
the wall, it "sticks" (is absorbed) , thus giving, as t-*- 00 , 
the distribution of win probabilities. (It is interesting 
to note that this relates the work of Clark [Ref. 8], 

Weiss [Ref. 35], Morse and Kimball [Ref. 22], and Snow 
[Ref. 27], who were concerned with the solution to the 
Chapman-Kolmogorov, to the work of Smith [Ref. 26], Brown 
[Ref. 7], etc., who were concerned with the win probabilities) . 
Note that in the present case, Y has a very large probability 
of winning. If the forces were moved closer to parity, more 
of the wave would be absorbed by the Y force boundary , 
giving X a larger probability of winning. 
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Plot of P(t,tn,n) for Fixed t 
FIGURE 1 
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P(t,m,n) 




Plot of P(t,m,n) for Fixed t 
FI CURE 2 
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P(t,m,n) 

I 




Plot of P(t,m,n) for Fixed t 
FIGURE 3 
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P(t,m,n) 




Plot of P(t,m,n) for Fixed t 
FICIJIIF 4 
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P(t t m,n) 




Plot of P(t,m,n) for Fixed t 
FIGU.1." 5 
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The above discussion has given a description of what 
occurs with the probability mass over time. Another des- 
cription .of the process is, it acts as a combination of 
diffusion and convective transport of probability mass. 

This concept enhances ones understanding of the nature of 
so-called diffusion approximations of Markov processes 
(i.e., state space assumed continuous). Additionally, it 
provides insights into temporal changes in variance in force 
levels (see Chapter III, Section C) . Thus, it is felt that 
computer graphics, if utilized more fully could provide 
easily obtainable insights into the dynamics of a combat 
model, that would be difficult to obtain otherwise. 

C. COMPARISONS 

Two types of comparisons have been made between stochas- 
tic and deterministic models; probability of winning, and 
the expected force levels against the time history of the 
force levels of the deterministic model. Two recent studies 
have been done comparing the probability of winning with 
the certainty of winning or losing in the deterministic 
model [Ref. 20 and Ref. 28], Both came up with basically 
the same conclusions. 

(1) As force levels increase, the probability 

of winning approaches the certain probability 
of winning or losing in the deterministic 
model . 

(2) Significant differences occur only when the 
forces are near parity. 
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Lee and Wannasilpa [Ref. 20] came to their conclusions for 
the square-law, linear-law, and mixed-law. Springall arrived 
at the same conclusions for a slightly more complicated model. 

Snow [Ref. 27] showed that the time history of the 
expected force levels differed from the time history of 
the deterministic force levels; that is m(t)-x(t) ^ 0. 

Clark [Ref. 8] defined this as bias . This seems to be an 
appropriate term and will be used here. Snow apparently 
first developed an explicit expression for this bias for 
the square- law; 

n 

dm — ° 

■sr = - an + a 1 P(t,0,n) ( 37 ) 

. ar n=l 

Clark [Ref. 8] empirically showed that bias exists for the 
linear-law. All researchers seem to have arrived at the 
same general conclusions : 

(1) As initial force levels increase, the bias 
decreases . 

(2) The longer a battle lasts, the greater 
the bias. 

Powers and Taylor [Ref. 24] and Hanna [Ref. 12] have 
compared differences in optimal time sequential fire distri- 
bution policies for deterministic and stochastic models: 
Hanna's investigation indicates that there is a real differ- 
ence between the two (at least for small numbers of total 
combatants) . 
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Thus, some work has been done to investigate the differ- 
ences between deterministic and stochastic models, especially 
in the areas of who is going to win and the average time 
histories of the force levels. However, the analysis done 
up to the present has been somewhat limited in scope and 
does not allow an analyst to say with any degree of cer- 
tainty, when the complex random process of combat can be 
adequately represented by a deterministic model. 
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III. COMPARISON OF DETERMINISTIC AND STOCHASTIC 



MODELS FOR "SQUARE LAW” ATTRITION PROCESS 

It would be most useful to be able to make comparisons 
between stochastic and deterministic models used today. 

This is not feasible within the scope of this paper for 
several reasons. 

1. It is not clear, in many cases, when the two models 
are equivalent. For example, one could try and 
compare DYNTACS-X and Bonder/I . U. A. , but before 
this comparison could be made, one would have to 
insure that the input parameters are equivalent. 

In the case of these two models, that is not a 
trivial task, as the structure of the two models 
is quite different, particularly in route selection. 
About the best that can be done is to construct the 
input parameters for each model and assume they 
are equivalent. 

2. Parametric analysis is not feasible in a complicated 
model. Complicated models, both stochastic and 
deterministic, are expensive to run. Because of 
this and the complex interactions among the variables 
of a model, parametric analysis is not only 
prohibitively expensive, but difficult to do. 

3. It is difficult to gain insights in a comparison 
of simple stochastic and deterministic models, and 
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essentially impossible to gain many useful insights 
in a comparison of complicated models. 

Lanches ter- type models may be thought of as resulting from 
an aggregation of the myriad of combat environmental varia- 
bles into constant attrition rate coefficients and time 
varying force levels and, moreover, as complementary complex 
system models. Such idealizations will be considered in 
this thesis with the knowledge they are idealizations, but 
that they may provide insights that can be used as points 
of departure for the investigation of more complex models 
of combat. Specifically, two models will be compared; the 
deterministic Lanchester square-law and the square-law 
stochastic attrition process. It is believed that trends 
found from this investigation may be carefully used to 
determine trends for a comparison of more complex models. 

As was shown in the literature survey, the square-law 
is a special case of the general differential equations 

gjr = “ F(t,x,y) (15) 

= - G ( t, x,y) (16) 



where 



F(t,x,y) = ay 



and 
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G(t,x,y) = bx 



giving 



dx 

dt 



= " a Y 



( 1 ) 



dt 



= - bx 



( 2 ) 



with initial conditions 



x ( t = 0) = x. 



(3) 



y(t = 0) = Y r 



(4) 



Also shown was the general form of the stochastic model: 



" '~ dV ' ~ = F ( t ,m+l , n) P ( t,m+l , n) + G(t,m,n+l)P (t,m,n+l) 



- (F(t,m,n) + G(t,m,n) ) P (t,m,n) (20) 



with boundary conditions 



dP ( t , m, ,n ) 



~l re E ' ~ bp '~ = F(t ' m bp +1,n bp )p(t '”bp +1 'V > 



+ Gft.i^p.n^+DPlt.ll^p.l^p+l) (38) 
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dP(t,m,n. ) 

^ 2- = F(t,m+l,n bp ) P (t,m+l,n bp ) 



+ G(t,m,n bp +1) P (t,m,n b +1) 



F(t ,m,n bp ) P(t / m,n bp ) 



m > m 



np 



dP ( t,m ,n) 

-rr^ = F(t,m +l,n) P (t,m +l,n) 

dt np np 



+ G(t,m ,n+l)P(t,m ,n+l) 
np np 



- G(t,m ,n)P(t,m ,n) 
np np 



n ' "b, 



and initial conditions 



P (0 ,m, n) = 



1 m = m and n 
o 

0 otherwise 



= n 



Substituting the square-law attrition rates gives the 
square-law stochastic attrition process: 



dP ( t , m , n ) _ 



dt 



= anP(t,m+l,n) + bmP ( t,m,n+l) - (an+bm)P(t 



and boundary conditions 



( 39 ) 



(40) 



(41) 



,m,n) 

( 19 ) 
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dt 



an bp P (t ' m bp + 1 ' n bp 1 + bm bp P(t ' m bp'V +1) 

(41) 



an^P ( t,m+l,n bp ) + bmP (t, 111 , 1 ^+ 1 ) 

" an bp P{t,m,n bp ) (42) 

anP (t , m ri p + 1 ^ n ) + bi^P ( t,^ ,n+l) 

- bln bp P(t ' m np' n) (43) 

and initial conditions 

t 

P(o,m,n) = < 

,0 

( 

It was noted that the square-law comparisons have been 
made previously. However, as far as the author can tell, 
no previous comparisons have been made contrasting the models 
for breakpoint rules other than annihilation. Taylor [Ref. 30] 
has pointed out (and Adkins [Ref. 1] examined in further 
detail) that a neglected area in the modeling of land combat 
is the modeling of battle termination. In actual battles, 
several different events might cause battle termination. 

If one of the opposing forces is annihilated, then the battle 
ends; however, this is a rare event essentially never observed 



m = m and n = n 
o o 

otherwise 



(44) 



dP (t,m np ,n) 
dt 



dP(t,m,n bp ) 

dt 
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in real combat, with a few exceptions (e.g., Iwo Jima, Alamo, 
etc.). If one of the opposing forces surrenders unconditionally, 
the battle will also terminate. The third and most common 
event which will result in battle termination is that of one 
of the opposing force breaking contact with the enemy and 
withdrawing from the battlefield [Ref. 13). There are other 
possible events which might result in battle termination, 
but they will not be considered. 

For the purposes of further discussion, a breakpoint is 
defined to be that state of a battle at which a unit con- 
siders itself no longer capable of performing its mission 
and as a result elects to break contact with the enemy and 
withdraw from the battlefield. Therefore, when a unit with- 
draws from the battlefield strictly to avoid further combat, 
the unit is considered to have reached its breakpoint and 
lost the engagement. 

In many current land combat models a unit's breakpoint 
is determined by the percentage of casualties sustained by 
that unit [Ref. 1] . Although available empirical evidence 
shows that a deterministic battle termination model does not 
appear reasonable, it is an adequate "first cut" [Ref. 13]. 

And a breakpoint determined by the percent of casualties 
sustained is more reasonable than one that assumes the battle 
goes to annihilation. For this reason, the investigation 
was done with deterministic breakpoints determined by the 
percent of casualties. 
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A. PROBABILITY OF WINNING 



Hypothesis 1-1 



i 



As the breakpoint force levels, ^ and n, are 
moved closer to the initial forcerevels ,m o and 
n , respectively; the difference in the 
probability of winning between the two models 
increases . 

Lee and Wannasilpa [Ref. 20] deduced that as the initial 
force levels are increased, the difference between the proba- 
bility of winning in the deterministic model and the stochas- 
tic model decreases. In fact it appeared that the difference 
disappears in the limit. However, the conclusion reached 
by Lee and Wannasilpa was based only on battles to annihila- 
tion. A look at the recursive relationship used in Appendix 
A to solve for the probability of winning is based on a 
summation from 1^+1 to m Q . Furthermore, the recursive 
relationships are from m^ to m Q . These seem to indicate 
that the differences in the probability of winning decrease 
not just as m Q increases, but as m^-m^p increases. The 
following figure seems to support this. Note, that for a 
deterministic model, the probability that Y wins, P^, has 

a discontinuous jump from one to zero as the ratio 
bm o 2 (l-(l-F x ) 2 ) 

5 tt — passes through one. (Note: F and F are 

an Q ^ ( 1- (1-F y ) j X 

the fractions of initial combatants each 

side is willing to lose.) For the square- law stochastic 

attrition process, however, Py moves continuously from one 

to zero reaching about 0.5 at a ratio of one. Note that as 



F^ and Fy decrease, m^ and n^ increase, so m 0 -m b 



and 
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P v for Varvinrt breakpoints 
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n -n, decrease. As they decrease, P„ for the stochastic 
o bp J Y 

model moves further and further from Py for the deterministic 
model. This was found to be true for all combinations of 
parameters considered. 

Hypothesis 1-2 

For fixed breakpoint casualty /initial force 
ratios (F and Fy) , the difference between 
the probability of winning for the deter- 
ministic model and the stochastic model 
decreases as the initial force levels increase. 

Lee and Wannasilpa [Ref. 20] deduced the above except 
their conclusion was restricted to battles to annihilation. 
The above hypothesis generalizes the conclusion to include 
any set of fixed F^ and Fy. Once hypothesis 1-1 along with 
Lee and Wannasilpa' s conclusion are accepted, hypothesis 
1-2 follows. The following figures support this. It is 
interesting to note that the plots of Py are similar for 
all fixed sets of F^ and Fy. The only difference is in the 
angles of the curves. As F^ and Fy increase, the result 
is similar to rotating the curves counter-clockwise. This 
causes an increase in the difference between Py for the two 
models . 

Hypothesis 1-3 

As the forces move away from parity, the 
difference between the probability of 
winning for the two models becomes negligible. 

This was the final conclusion reached by Lee and Wanna- 
silpa. Although their conclusion was for battles to 
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P(Y wins) 
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annihilation, the present findings do not alter the conclu- 
sion for other fixed force level breakpoints. The only 
change is in how far from parity the forces must be to make 
the difference negligible. 

The general consequence of this hypothesis is that unless 
the forces are near parity, there is little difference in 
the probability of winning between the two models. From 
this viewpoint, then, if one side is going to win decisively, 
there is little reason to use the more expensive stochastic 
model. Indeed, this has some intuitive appeal: The more 

decisively a side wins the battle, the less influence any 
uncertainties will have on the outcome. On the other hand, 
if the forces are near parity, the uncertainties in the 
evolution of the dynamics of combat may cause the expected 
loser to win. 

If the previous hypotheses are combined with this hypothe- 
sis, some interesting trends can be observed. It has been 
concluded that, with increasing m^-m^p and n^n^ , the 
differences in the probability of winning between the two 
models decreases. When combined with hypothesis 1-3, this 
says that, with casualties on each side on the order of 20 
or more, and with forces not near parity, there is no 
significant difference between the probability of winning 
for each model. Even if the forces are near parity, if the 
force levels are very large, and both sides are willing to 
take 20 or more casualties ■, the difference in model win 
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probabilities is negligible. In these cases, if only the 

win probabilities were considered as measurement criteria 

between the two models, one would use the less expensive 

deterministic model as an acceptable substitute for the 

more realistic stochastic model. 

One must be cautious of this conclusion. It would be 

inaccurate to make this conclusion unless the only objective 

of the model is to determine which side will win. In most 

cases, there are other insights analysts try to gain from 

the model, such as how many casualties does each side take, 

at what points in the battle are the most casualties being 

taken, how long will the battle last, etc. 

Nevertheless, the conclusions of this section do give 

indications of when the two models will not give the same 

winners. In these cases, no further evaluation of the models 

is needed, because it is obvious the deterministic model is 

not an acceptable substitute. As Springall [Ref. 28] said; 

Probably the single most important 
comparison is between the two predic- 
tions of which side is going to be 
victorious. If the two methods cannot 
agree on this, there is little hope 
for agreement on the subsidiary 
attributes . 

B. TIME HISTORY OF EXPECTED FORCE LEVELS 

If a determination has been made that the models will 
predict the same winner, further analysis may be made to 
see how well the models compare. Clark [Ref. 8] suggests 
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that a useful comparison is the time history of the expected 

number of survivors. In many models, one of the most utilized 

variables (performance and/or proxy variables) in measures 

of effectiveness for systems being evaluated by the models 

is the number of survivors. It has been stated previously 

that a bias exists in the square-law and the linear- law. 

Lanchester himself suspected as much; he recognized that 

his differential equations were approximations to the casualty 

rates likely to be experienced in an actual battle. Quoting 

from Lanchester [Ref. 19]; 

Since the forces actually consist of a 
finite number of finite units (instead of 
an infinite number of infinitesimal units) 
the end of the curve must show discontinuity, 
and break off abruptly when the last man is 
reached; the law based on averages evidently 
does not hold rigidly when the numbers become 
small. 

Lanchester stated that his differential equations are based 
on averages implying an underlying stochastic process. He 
also suggests that his differential equations may be good 
approximations only so long as the force sizes are large. 
Although his reasoning is heuristic, Snow [Ref. 27] verified 
that the square-law is biased for small numbers. He showed 
that 



dm 

dt 



-an + 



a 



S 

n=l 



nP ( t , 0 , n) 



m 

dn ^ 

= -bm + b Z mP(t,m,o) 

m=l 



(37) 



(45) 



53 



The conclusion is that the solutions to Lanchester's equa- 
tions for the square- lav; are good approximations to expected 
survivors so long as 

n 

o 

a Z nP(t,0,n) ~ 0 (46) 

n=l 



b E mP(t,m,0) ~ 0 (47) 

m=l 

Intuitively this seems sound. (46) and (47) express the 
boundary effects at time t. At early times in the battle, 
the expected force levels are near the initial force levels, 
and there is very little of the probability mass at the 
boundary. Thus the boundary effects are small. As the 
battle nears termination, though, much of the probability 
mass is near the boundary and (46) and (47) will be signi- 
ficant enough to effect the process. As there are no boun- 
dary effects in the deterministic process, there will be 
a difference in force level trajectories of the two processes. 

Clark provided a derivation of (37) and (45) . Additionally, 
he provided the following argument showing that the linear- 
law, which at first appears to be unbiased, is not. Through 
a simple derivation, it can be shown for the linear-lav;: 

Trt = - aE(MN:t) (48) 
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(49) 



gir = - bE(MN:t) 

There are no bias terms in this expression (as opposed to 
(37) and (45)) since the random variable MN has a value of 
zero whenever M = 0 or N = 0. Consequently, the assumed 
differential equations for the linear-law appear unbiased. 

However, even though the analysis of the solutions to 
these equations show that the state solution is unbiased, 
the time solution interjects a source of bias. This is 
easily seen as a result of elementary probability theory: 

E (MN) = E (M) • E (N) if an only if M and N 

are independent . 

In the linear-law, they are obviously dependent on each 
other. The rate of change of each force level is dependent 
on the product of the two force levels, causing the force 
levels to be dependent. 

Clark [Ref. 8] also gave numerical examples to show the 
amount of bias resulting from the difference in models. His 
results showed that there are significant amounts of bias 
resulting in both the linear-law and square- law cases. 

There are two things that should be noted about his results : 
First, his results are for small numbers of combatants; and 
second, his results are only for battles to annihilation. 

It was not shown in his dissertation that there is a 
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significant bias with large numbers of combatants or with 
breakpoints other than annihilation. 

The following hypotheses are based on extension of his 
conclusions to include larger forces and fixed force level 
breakpoints. Appendix D is a derivation of the square-law 
stochastic attrition process showing that the bias noted 
by Snow exists with any set of fixed breakpoint force levels 
Unfortunately, the form of the bias is not as clear as that 
for annihilation, since there are two additional terms. The 
use of the graphic aids noted earlier gives some indication 
that non-zero breakpoint force levels do not change the 
general progress of the battle. A look at the following 
figures will show that the time state probabilities progress 
in the same manner as in the case of battles to annihilation 
Figures 12-15 are plots of the time state probabilities of 
the same model shown earlier except the breakpoint force 
levels are eight instead of zero (annihilation) . Figures 
16-19 are the same except the breakpoint force levels are 
sixteen. The differences are that the "wave" is absorbed 
by the boundary at (m^^n) and (rt^n^) instead of at (0,n) 
and (m,0) . The question now is. What is the effect on 
the bias? 

A first analysis of (101) and (102) might give the 
indication that the bias is always positive (i.e., 
m(t)-x(t)> 0 for all t > 0) . However, this is not true, 
as can be seen in the following argument. 



56 



P(t,n,n) 




Plot of P(t,m,n) for Fixed t 
FICIJPu" 12 
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P(t,m,n) 




Plot of P(t,m,n) for Fixed t 
FIGURE 13 
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P(t ,m,n) 




Plot of P(t,m,n) for Fixed t 
FIGURE 14 
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P(t,m,n) 




Plot of P(t,m,n) for Fixed t 
FICURE 15 
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Plot of P(t,m,n) for Fixed t 



FIGURE 16 



P(t,m,n) 




Plot of P(t,m,n) for Fixed t 
FIG M. In 17 
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P(t,m,n) 




\ 




Plot of P(t,m,n) for Fixed t 
FICUSS 18 
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P(t,m,n) 




Plot of PCt^m^n) for Fixed t 
FICUSE 19 
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Define 



A x = m(t) - x ( t) 



Ay = n ( t) - y ( t) 



then 



A v ( t-0 ) = m(t=0) -x ( t=0) = - an + an = 0 

A O O 



Ay(t=0) = n(t=0) -y (t=0) = - bm Q + bm Q = 0 



Let 



s Y (t) 



n 

a[ £ nP(t ' m np' n) + ”bp I V P<t ' , ”np +1 ' n bp ) 

n n bp 

m 

+ v m £, + 1 p(t ' m ' n b P n 

m=m +1 c 

np 



m 



s x (t) 



b[ I mP(t,m,n > + n h p n bp P * * ' m : ,n +1) 
m=m * r r l 



np 



n 



+ m l P (t ,m. ,n) ] 

np . -i bp 

^ n=n, +1 L 



bp 



If 



(50) 

(51) 

(52) 

(53) 

(54) 

(55) 
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ay 



( 1 ) 



and 



then 



dx 

dt 



= - an + S y (t) 



dm _ dx 
dt dt 



- an + S (t) + ay 



d (m-x) 
dt 



- a(n-y) + S y (t) 



Substituting (50) 



dA 



X 



dt 



“ aA y + S y (t) 



and equivalently 



dA y 

"dt" 



- bA x + S y (t) 



with initial conditions 



A x (t = 0) = 0 



( 101 ) 



(56) 



(57) 



(52) 
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Ay ( t = 0) = 0 



(53) 



Define 



A = 



*x\ 

V 



and U = 



/s v (t>' 



s x (t) 



then 



d~ (° a 

dt = " , ^ 

b o 



:a+ u(t) 



and the solution of (58) is given by 



where 



A ( t) = e“ At A(t=0) + / e" A(t “ T) U(T) d T 



A = 



0 a 



i b °j 



Proof : 



d£ 

dt 



— A +- 

-Ae A(t=0) + U(t) 



A / e A(t T) U( T )d T 



= -AA (t) + U(t) 



(58) 



(59) 



AA(t) 



67 



from Sylvester's 



Theorem 



1 



e 



-AT 



cosh ( "V ab t) 
^-'yb/a sinM'yab t) 



- ~\J a/b s inh ( \ ab t) 
cosh ( ■y/ ab t) 



Hence 



A x ( t ) 



A y (t) 



t 

/ S (t) coshCy ab (t-x)) 

0 



Sx( T )“\ a/ b sinh(^ab (t-x)) ax 

(60) 



t 

f S (x)cosh(yab (t-x)) 
0 X 



- Sy ( x ) ~\J b/a s inh Cyab (t-x))dx 

(61) 



If one 

(m 

np 



now takes the case where one force is annihilated 
n b = 0 ) • then S^t) and S^(t) simplify to 



n 

o 

S (t) = E nP (t ,0 ,n) 
1 n=0 



S v (t) = E mP(t,m,0) 
X m=0 



(62) 



(63) 



which are simply the average number of Y survivors when X 
is annihilated at time t and the average number of X sur- 
vivors when Y is annihilated at time t respectively. It 



Sylvester's Theorem: For a polynomial function or 

convergent power series of a matrix, when the matrix has 
distinct eigenvalues. A., 

n n^ A-A.I 

F(A) = E F(A.) E (-r 3 — ) where A is an nxn matrix. 

i=l j-1 i j 

j^i 
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seems intuitively appealing that if X wins decisively, 

S Y (t) ->-0. If this is true, then (60) will become negative 

and (61) will become positive. This seems to be confirmed 

by numerical results obtained. 

Another interesting case to examine is the case of 

complete parity (a=b, m o =n Q , m np =n bp ) where m np =n bp =0 • 

In this case, S (t) = S (t) = S(t) — J 

x A 



A v (t) = / S ( t) [cosh ( V ab ( t — t ) ) - sinh ( y ab (t-x))]dx 
Y 0 

(65) 

By inspection, A^(t) > 0 and A y (t) >0 for t > 0, and both 

biases are positive. This is reinforced by numerical 
results . 

The general conclusion reached from this analysis is 
that the biases are both positive unless one side is going 
to win very decisively. In this case, the bias on the 
winner's side will be negative, while the loser's is positive. 

The remainder of this section will be a statement and 
discussion of hypotheses concerning bias. 

Hypothesis 2-1 

Given fixed initial force levels and fixed 
attrition coefficients, as the breakpoint 
force levels increase, the numerical bias 
decreases. But as a percentage of casualties 
of the deterministic model, the bias increases. 




A (t) = / S ( x) [cosh ( "W ab (t-x)) - sinh ( \/ ab (t-x))]dx 

0 

(64) 



t 



t 
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This hypothesis has heuristic appeal when the bias 
equations are analyzed. The bias terms are (see Appendix D) : 



n 



S Y (t) 



a '_L nP(t ' m np' n) + IVV (t ' In np +1 ' n bp ) 



n=n 



bp 



m 



p bp 



+ n Z P ( t ,m,n ) ] 

bp m=m +1 bp 

np 



( 54 ) 



m 



s x (t) 



ME + ra n p n bp P<t ' m bp' n b p + l> 

m “bp 



n 



+ nr Z P (t ,m ,n) ] 

tp n=n, +1 np 

bp 



( 55 ) 



As n^ and increase, the number of terms in each summation 

decreases. Although the second and third terms are multi- 
plied by n^p (which is increasing) , it would seem that this 
may be more than offset by the decrease in number of terms. 

Even if this is true, as a percentage, the bias will 
increase unless the actual bias decreases as fast as ia^-m^ . 
Although parametric analysis cannot logically prove the truth 
of this statement, it shows that the statement is certainly 
plausible . 

As all parameters except m^p and n^p are fixed, a 

numerical parametric analysis can be performed. A range 

of fixed parameters was chosen; small and medium values of 

m and n , and several combinations of a and b. In all 
o o 

cases the results generally agreed. The following figures 
illustrate the results. They plot A % and A v % as defined in 

A X 
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Time History Hhowinc* Hias 
FIGURE ?0 
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% of Survivors on m Side of Total Casualties for Deterministic flattie 




a=0*008 
b=0.004 
M =40 
M? =8 
^40 

“bp* 8 




Percent of Ti'ne to Completion of Deterministic 
Battle 



Time History Showing Hias 
FIGURE 21 
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Battle 



Time History Showing Rias 
FIGURE 22 
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Battle 



Time History Showing Mias 
FISH RE 24 
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Time Vistory Showing ias 
FIGURE 25 
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Battle 



Time History Showing Rias 
FIGURE 26 
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Battle 



Time History Showing Mas 
FIGURE 27 
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Battle 



Time history Showing ^ias 



FIGURE 2°s 
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Battle 



Time History Showing P.ias 
FIGURE 29 
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A Survivors 




Mbp-^bp 




M bp= N bp 
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Appendix C as functions of scaled time as defined in the 
same appendix. The important things to notice are when the 
lines separate and how large the final separation is for 
the various input parameters listed on each figure. The 
figures are summarized in Figures 30 and 31. 

In our examples, X is taken always to be the loser and 
Y to be the winner. As the breakpoints are increased, the 
difference in X survivors decreases from slightly over three 
to slightly over one. This seems to confirm the first part 
of the hypothesis. Additionally, the rate of change of the 
bias is less than the rate of change of the breakpoints, so 
the percentage of bias increases, as indicated in Figure 31. 
This seems to confirm the second part of the hypothesis. 
However, the bias for the winning side does not follow the 
hypothesis (Y wins very decisively in this case) . Indeed, 
it starts negative, becomes positive, and then decreases. 

But two things must be noted before the hypothesis is re- 
jected: First, the slope of the Ay curve is small; and 

second, it eventually decreases. This leads to the Ay% 
curve in Figure 31, where Ay% is increasing. 

It is argued in Appendix C that A% is the more important 
measure of the bias. If this is accepted, then the second 
half of the hypothesis is more important and, as indicated 
above, it is true. All of the results in Table 2 reinforce 
this . 

One final way of stating this conclusion is; the greater 
the differences m 0 -rn bp and n 0 rn bp' t ^ ie ^ ess significant the 
biases . 



82 



a 


b 


m 

o 


n 

o 


Uk. 


3* 




V 




V 


0.008 


0.004 


40 


40 


0 


0 


3.22 


8.05 


- 0.23 


- 1.96 


0.008 


0.004 


40 


40 


8 


8 


2.62 


8.19 


0.18 


1.61 


0.008 


0.004 


40 


40 


16 


16 


2.11 


8.79 


0.42 


4.40 


0.008 


0.004 


40 


40 


24 


24 


1.64 


10.25 


0.52 


7.41 


0.008 


0.004 


40 


40 


32 


32 


1.15 


14.38 


0.47 


12.43 


0.008 


0.004 


24 


24 


0 


0 


2.48 


10.33 


- 0.22 


- 3.13 


0.008 


0.004 


24 


24 


4 


4 


2.10 


10.50 


0.06 


0.88 


0.008 


0.004 


24 


24 


9 


9 


1.68 


11.20 


0.29 


4.93 


0.008 


0.004 


24 


24 


14 


14 


1.31 


13.10 


0.38 


8.74 


0.008 


0.004 


24 


24 


19 


19 


0.94 


18.80 


0.38 


16.10 


0.004 


0.0015 


40 


40 


16 


16 


2.06 


8.58 


0.28 


4.06 


0.004 


0.0015 


40 


40 


24 


24 


1.62 


10.13 


0.36 


7.02 


0.004 


0.0015 


40 


40 


32 


32 


1.12 


14.00 


0.33 


11.79 


0.004 


0.0015 


24 


24 


0 


0 


2.28 


9.50 


0.13 


- 2.58 


0.004 


0.0015 


24 


24 


4 


4 


1.98 


9.90 


0.04 


0.82 


0.004 


0.0015 


24 


24 


9 


9 


1.63 


10.87 


0.19 


4.48 


0.004 


0.0015 


24 


24 


14 


14 


1.28 


12.80 


0.27 


8.49 


0.004 


0.0015 


24 


24 


19 


19 


0.90 


18.00 


0.25 


14.37 



TABLE II. Bias for Different Breakpoints 
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Hypothesis 2-2 

For fixed and F^, and fixed attrition 
rate coef f icients , A the larger the initial 
forces, the greater the numerical bias, 
but the smaller the percentage bias. 

A generally accepted conclusion for Lanchester-type 
models is a type of "law of large numbers". For a battle 
to annihilation, the differences between a deterministic 
model and the equivalent stochastic model "disappear" as 
the initial force levels increase. Hypothesis 2-2 is a more 
general statement of this not for battles to annihilation, 
but to fixed casualty/initial force level breakpoints. 
Additionally, it is concerned only with the differences in 
force level histories. Hypothesis 1-2 was a discussion of 
the "law of large numbers" for the probability of winning. 

There is an intuitive appeal for the hypothesis. The 
time history of the force levels for a deterministic model 
is something of an "average" time history of the stochastic 
model. As the number of combatants involved increases, the 
greater the aggregation, and the closer this average is to 
the mean force level history of the stochastic model. 

All but two parameters are held constant, so a numerical 
parametric analysis can easily be done. Figures 31 to 35 
show some of the results of the variation in parameters. 
Figures 37 and 38 summarize them and show that, for these 
fixed parameters, the hypothesis is true. 

Graphs of the results in Table III would give similar 
results and further reinforce the hypothesis. 
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a 


b 


fx 


fy 


m 

o 


n 

o 




V 




V. 


0.008 


0.004 


1.0 


1.0 


12 


12 


1.75 


14.58 


-0.20 


- 5, 


0.008 


0.004 


1.0 


1.0 


24 


24 


2.48 


10.33 


-0.22 


- 3, 


0.008 


0.004 


1.0 


1.0 


40 


40 


3.22 


8.05 


-0.23 


- 1 , 


0.008 


0.004 


0.8 


0.8 


24 


24 


2.10 


10.50 


0.06 


0 , 


0.008 


0.004 


0.8 


0.8 


40 


40 


2.62 


8.19 


0.18 


1 . 


0.008 


0.004 


0.6 


0.6 


24 


24 


1.68 


11.20 


0.29 


4. 


0.008 


0.004 


0.6 


0.6 


40 


40 


2.11 


8.79 


0.42 


4. 


0.008 


0.004 


0.4 


0.4 


24 


24 


1.31 


13.10 


0.38 


8. 


0.008 


0.004 


0.4 


0.4 


40 


40 


1.64 


10.25 


0.52 


7. 


0.008 


0.004 


0.2 


0.2 


24 


24 


0.94 


18.80 


0.38 


16, 


0.008 


0.004 


0.2 


0.2 


40 


40 


1.15 


14.38 


0.47 


12. 


0.004 


0.0015 


1.0 


1.0 


6 


6 


1.13 


18.83 


-0.10 


- 7. 


0.004 


0.0015 


1.0 


1.0 


8 


8 


1.31 


16.38 


-0.11 


- 6. 


0.004 


0.0015 


1.0 


1.0 


12 


12 


1.61 


13.42 


-0.12 


- 4, 


0.004 


0.0015 


1.0 


1.0 


24 


24 


2.28 


9.50 


-0.13 


- 2. 


0.004 


0.0015 


0.6 


0.6 


24 


24 


1.63 


10.87 


0.19 


4. 


0.004 


0.0015 


0.6 


0.6 


40 


40 


2.06 


8.58 


0.28 


4, 


0.004 


0.0015 


0.4 


0.4 


24 


24 


1.31 


13.10 


0.38 


8. 


0.004 


0.0015 


0.4 


0.4 


40 


40 


1.64 


10.25 


0.52 


7. 


0.004 


0.0015 


0.2 


0.2 


24 


24 


0.94 


18.80 


0.38 


16. 


0.004 


0.0015 


0.2 


0.2 


40 


40 


1.15 


14.38 


0.47 


12. 



TABLE III. Bias for Different Initial Force Levels 
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The general result of this hypothesis is that, if the 
forces are of any appreciable size (greater than 20) , with 
breakpoint force levels such that in^-m^p and n^n^ both 
exceed 20, the deterministic model is probably an adequate 
representation of the complex random process of combat. 

Hypothesis 2-3 

If any parameter is varied to bring the 
forces closer to parity, the time of the 
battle increases and the bias increases. 

This hypothesis is a result of an analysis of the solu- 
tion for A^(t) and Ay(t). As was shown previously, 

t 

Ay ( t) = / [S (x) cosh( Yab(t-x) ) -Sy(x) \ a/bsinh(\ ab(t-x) ]dx 

0 

(60) 

t 

Ay ( t) = f [S x ( x) cosh ( -yaF( t-x) ) -S y ( x) a/bsinh ( Y ab (t-x) ] dx 

(61) 

It appears that, as t increases, the integral will increase 
in absolute value. As only one parameter is varied, numeri- 
cal parametric analysis was easy. All results computed 
indicated the hypothesis is true. Figure 40 shows typical 
results. It is interesting to note that the hypothesis 
holds for the loser. However, it is not clear what happens 
for the winner, but it would appear that it is a continuous 
function. If it is, then at some point, A^ reaches a 
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maximum negative value, and then increases to a positive 
value at parity. It then appears the hypothesis only holds 
for the loser. The results in Table IV show the same 
trends for the winner and loser. 

Hypothesis 2-4 

If a battle were to be terminated at a 
time t<0.5tf, where t^ is the completion 
time of the deterministic model, there is 
no significant difference between the 
results of the two models. 

Decisions to terminate battles are not always made 
because too many casualties have been taken. A good example 
of this would be the rear guard action of a division involved 
in a retrograde movement. The rear guard delays the enemy 
force, in many cases, long enough to allow the main force 
to retrograde and set up a defensive position, or at least 
allow a fresh rear guard to prepare for a further delaying 
action. In either case the rear guard commander will be 
tasked to delay the enemy for a fixed amount of time. When 
that time limit is reached, he will disengage, regardless 
of whether he feels he could engage the enemy longer. In 
fact, he may have taken very few casualties. This would be 
a situation in which hypothesis 2-4 would be applicable. 

An evaluation of A^(t) and A^(t) will show that they 
are generally increasing (in absolute terms) , with increasing 
time. A study of various combinations of parameters has 
shown this to be the case. The question then is, at what 
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a 


b 


o 


0.008 


0.004 


24 


0.008 


0.004 


24 


0.008 


0.004 


24 


0.008 


0.004 


24 


0.008 


0.004 


24 


0.004 


0.001 


6 


0.004 


0.001 


12 


0.004 


0.001 


6 


0.004 


0.0015 


6 


0.004 


0.002 


6 


0.004 


0.004 


6 


0.004 


0.004 


6 


0.002 


0.004 


6 


0.0015 


0.004 


6 


0.001 


0.004 


6 



?k 


n 

o 




fbc 


19 


24 


0 


0.88 


19 


24 


0 


1.30 


9 


24 


0 


1.68 


4 


24 


0 


2.10 


0 


24 


0 


2.48 


0 


6 


0 


1.06 


0 


6 


0 


3.22 


0 


6 


0 


1.06 


0 


6 


0 


1.13 


0 


6 


0 


1.23 


0 


6 


0 


1.93 


0 


12 


0 


1.06 


0 


12 


0 


1.18 


0 


12 


0 


1.34 


0 


12 


0 


1.97 



bm 2 (1-(1-F ) 2 ) 
o x 

A Y t f an o 2 (l-(l-F Y ) 2 ) 



0.35 


27 


0.18 


0.37 


58 


0.32 


0.29 


91 


0.42 


0.06 


126 


0.48 


-0.22 


156 


0.50 


-0.07 


275 


0.25 


1.93 


493 


1.00 


-0.07 


275 


0.25 


-0.10 


291 


0.375 


-0.14 


312 


0.50 


1.93 


CO 


1.00 


-0.11 


137 


0.25 


-0.32 


312 


0.50 


-0.47 


493 


0.75 


3.22 


CO 


1.00 



TABLE TV. Bias for Different Attrition Rate Coefficients 
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time into the battle do A (t) and A (t) become significant? 
No method was found to determine this, but Figures 41 to 69 
are typical of all results found. 

As might be expected, the time at which the two time 

histories begin to diverge is strongly related to the final 

differences A^%(t^) and Ay%(t^) . As has been noted, as 

m -m. and n -n, decrease, and the force levels are moved 
o bp o bp 

closer to parity, A x %(t^) and A^% (t^) increase. The worst 
examples found did not diverge before 0.5 t^, with one 
excaption. The one exception to this hypothesis is when 
the forces are at or very near parity. In this case, the 
hypothesis does not hold because, with the forces near or 
at parity, the rate at which casualties on both sides occur 
approaches zero as time increases. This gives, in a sense, 
an exaggerated time length of the battle and thereby distorts 
the relative amount of time A^% and A^% are appreciable. 

C. VARIANCE 

It is possible for the two models to agree on which side 
is going to win, and on the time histories of the average 
force levels, but the deterministic model may not adequately 
represent the complex random process of combat. This is 
because the outcome of the stochastic model is determined 
by the outcome of many uncertain events. The probability 
distribution function related to these uncertainties may 
cause the variance of the force levels to be considerable . 
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If the system being evaluated includes measures of effec- 
tiveness that are functions of the distribution of the final 
force levels, then in the case of large relative variance, 
the deterministic model is not an adequate representation 
of the complex random process of combat. 

Clark [Ref. 8] analyzed the variance of the square-law 
stochastic attrition process and hypothesized there are 
two types of time histories of the variance. The first is 
where the variance grows and converges asymptotically to a 
final level, while the second is where it grows to a peak 
level and then decreases asymptotically to a final level 
(see Figures 70 and 71) . 

The qualitative pictures of the time state probabili- 
ties in Chapter II provide some useful insights into why 
the variance asymptotically approaches a value for every 
set of input parameters. As was noted, the change in the 
time state probability is a combination of diffusion and 
convective transport of the probability mass. The diffusion 
provides the variance in the force level. As the probability 
mass "hits" the boundary it remains there (is absorbed) ; 
for the probability mass absorbed, the diffusion effect 
has also ended and the variance of that probability mass 
is fixed. As more is absorbed, the less relative impact the 
remaining probability mass will have on the variance. Even- 
tually, as time grows, the remaining probability will have 
no effect, thus giving the asymptotic characteristic. 
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However, because of the convective transport of probability 
mass, the remaining probability mass will not only continue 
to be diffused, it will, in general, arrive at the boundary 
closer to the breakpoint values (e.g., for the X force, the 
center of mass of the small probability mass, AP^ arriving 
in the interval (t,t+At^), will be at m survivors, but 
AP 2 arriving in the interval ( t+At^ , t+At 2 ) will be centered 
at m-k survivors). Thus, as time increases, the convective 
transport of the probability mass causes the average of the 
probability mass not absorbed to approach the breakpoint 
force level. The breakpoint force level acts as a boundary, 
and the probability mass cannot be diffused past that point. 
In fact, the probability that would have been diffused past 
the breakpoint "piles up" at it, and the variance is not 
quite as great as it might be. In fact, if the probability 
mass is transported fast enough toward the breakpoint, the 
variance of the remaining probability mass may decrease in 
spite of the diffusion and might cause the overall variance 
to decrease. 

For illustration, take Figure 3 and collapse the 
probability mass to the X (m survivor) axis. The result 
is : 
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This gives a qualitative picutre of P(t,77.5,m) . The vari 
ance is calculated from the sum of the squared deviations 
from the mean. In this case, it appears the variance will be 
"large” (in fact, it is on the order of magnitude of the mean) 
If we now take Figure 5 and collapse the probability 
mass to the X (m survivor) axis, the result is: 



P(t=155.5,tn) 




C'M.I.AP^'O <’RUtVr.ILI nn Y MASS 



T .'i;*ure 73 



Obviously the variance is much small than it was at 
t = 77.5. 

Note that X is the decisive loser in this battle; the 
convective transport has moved the probability mass toward 
the breakpoint fast enough to offset the diffusion effect, 
and the variance has decreased. On the other hand, if the 
same graphs were made for the Y force (the decisive winner) 
it would be obvious that the probability mass is not being 
transported toward the Y breakpoint very rapidly, and the 
variance continues to increase. 

Thus it seems that how the variance progresses is depen- 
dent on how rapidly the side is being attrited, relative to 
the other side. If a side wins, the variance will grow 
asymptotically, whereas if a side loses, the variance will 
peak and then decrease asymptotically. The results of all 
combinations of parameters in this study indicate this type 
of performance. 

Furthermore, the results in Table V indicated two 
general trends . 

Hypothesis 3-1 

As m -m, increases, the variance increases 
in numerical value, but decreases in relative 
measure to the expected casualties. 

Hypothesis 3-2 

As the forces approach parity, the variance 
increases . 
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a 


m 

o 




b 


n 

o 
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"V V ( Y ) 
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0 
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40 


0 
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4.67 


0.008 


40 


8 


0.004 


40 


8 


3.66 


4.18 


0.008 


40 


16 


0.004 


40 


16 


2.97 


3.61 


0.008 


40 


24 


0.004 


40 


24 


2.32 


2.94 


0.008 


40 


32 


0.004 


40 


32 


1.71 


2.04 


0.008 


24 


0 


0.004 


24 


0 


3.45 


3.59 


0.008 


24 


4 


0.004 


24 


4 


2.90 


3.27 


0.008 


24 


9 


0.004 


24 


9 


2.34 


2.81 


0.008 


24 


14 


0.004 


24 


14 


1.85 


2.26 


0.008 


24 


19 


0.004 


24 


19 


1.35 


1.56 


0.008 


12 


0 


0.004 


12 


0 


2.39 


2.49 


0.008 


24 


0 


0.004 


24 


4 


2.90 


3.27 


0.008 


24 


0 


0.004 


24 


9 


2.34 


2.82 


0.008 


24 


0 


0.004 


24 


14 


1.82 


2.29 


0.008 


24 


0 


0.004 


24 


19 


1.33 


1.65 


0.004 


40 


16 


0.0015. 


40 


16 


2.89 


2.94 


0.004 


40 


24 


0.0015 


40 


24 


2.29 


2.46 


0.004 


40 


32 


0.0015 


40 


32 


1.68 


1.81 


0.004 


24 


0 


0.0015 


24 


0 


3.17 


2.79 


0.004 


24 


4 


0.0015 


24 


4 


2.73 


2.60 


0.004 


24 


9 


0.0015 


24 


9 


2.28 


2.29 


0.004 


24 


14 


0.0015 


24 


14 


1.77 


1.92 


0.004 


24 


19 


0.0015 


24 


19 


1.28 


1.40 


0.004 


12 


0 


0.0015 


12 


0 


2.20 


1.96 



TABLE V. Standard Deviation of Force Levels 
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The second trend is in general agreement with previous 
conclusions in the respect that, unless the forces are near 
parity, the models are not significantly different. The 
first trend is in agreement with Hypothesis 2-1, so that 
the conclusions reached for it seem to apply also for the 
variance. Thus, even though the variance is a different 
measure for comparison of the two models, it seems to agree 
with previous conclusions. 
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IV. SUMMARY 



The intent of this analysis was to consider the question, 
can the complex random process of combat be adequately 
represented by a deterministic model. The models chosen 
for comparison were the deterministic and stochastic attri- 
tion versions of the Lanchester square- law attrition pro- 
cess. Three aspects of these processes have been compared 
and analyzed. 



Deterministic Model Stochastic Model 



A. Fixed winner and loser 

B. Time history of the 
force levels 

C. No variance of the 
force levels 



Probability of winning 

Time history of the 
expected force levels 

Variance of the force 
levels 



Conclusions in the form of hypotheses were discussed for 
each of the above comparisons: 



A. PROBABILITY OF WINNING 
Hypothesis 1-1 

As the breakpoint force levels, m, and 
n, are moved closer to the initial 
force levels, m and n Q , respectively, 
the difference in the probability of 
winning between the models increases. 



Hypothesis 1-2 

For fixed breakpoint casualty/initial force 
ratios (F„ and F ) , the difference between 
the probability or winning for the deter- 
ministic model and the stochastic model 
decreases as the initial force levels increase. 
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Hypothesis 1-3 



As the forces move away from parity, the 
difference between the probability of 
winning for the two models becomes 
negligible . 



B. TIME HISTORY OF EXPECTED FORCE LEVELS 



Hypothesis 2-1 

Given fixed initial force levels and fixed 
attrition coefficients, as the breakpoint 
force levels increase, the numerical bias 
decreases. But as a percentage of casual- 
ties of the deterministic model, the bias 
increases . 



Hypothesis 2-2 

For fixed F^ and F y , and fixed attrition 
rate coefficients, the larger the initial 
forces, the greater the numerical bias, 
but the smaller the percentage bias. 



Hypothesis 2-3 

If any parameter is varied to bring the 
forces closer to parity, the time of the 
battle increases and the bias increases. 



Hypothesis 2-4 

If a battle were to be terminated at a time 
t £ 0.5 t^, where t^ is the completion time 
of the deterministic model, there is no 
significant difference between the results 
of the two models. 



C. VARIANCE OF FORCE LEVELS 



Hypothesis 3-1 

As m -itl increases, the variance increases 
in numerical value, but decreases in relative 
measure to the expected casualties. 
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Hypothesis 3-2 

As the forces approach parity, the variance 
increases . 

The hypotheses were formulated and evaluated through 
analyses of combinations of analytical solutions, numerical 
solutions, and graphical representations. 
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V. DISCUSSION 



The hypotheses formulated in this thesis provide a basis 
for the extension of several conclusions about the adequacy 
of the deterministic model. When one side is going to win 
decisively in the deterministic model, the models do not 
differ significantly and the deterministic model is to be 
preferred. How decisively one side must win to make the 
difference insignificant is a function of two sets of 
parameters; the initial force levels and the breakpoint 
force levels. As the initial force levels increase, the 
relative differences between models decrease as long as 
iu^-m^p and n 0 -n bp do not increase. Thus, with relatively 
large initial force levels (greater than 20) and breakpoint 
force levels such that there will be relatively large losses 
allowed (also greater than 20) , the models are not signifi- 
cantly different. 

It must be recognized that these conclusions only apply 
to the very simple and idealistic Lanchester square- law 
attrition process and the square-law stochastic attrition 
process. As is always the case, the insights into combat 
dynamics obtained from these models are no more valid than 
the models themselves. However, intuition indicates the 
conclusions may be generalized. It is intuitively appealing 
to say that, no matter what the model or degree of complexity, 
if the deterministic model predicts that one side will win 
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decisively, the stochastic model will not disagree 
significantly . 

Nevertheless, it is felt that the hypotheses may 
provide a point of departure for the comparison of more 
complex models of combat. 
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VI. CONCLUSIONS 



For the idealized model (i.e., the so-called Lanchester 
"square-law" attrition process) , it is concluded that the 
complex random process of combat can be adequately repre- 
sented by a deterministic model if the following conditions 
are met: 

1. Each side starts with at least 20 combatants. 

2. Each side is willing to take at least 20 casualties. 

3. The forces are not near parity. 

4. If the forces are near parity, but each side 
initially has in excess of 40 combatants and 
is willing to take in excess of 20 casualties. 

Consequently, it seems plausible that similar results would 
be obtained if one were to compare analogous deterministic 
and stochastic "real world" combat models (i.e., large- 
scale complex system models of military forces) . In other 
words, it seems reasonable to expect that if the opposing 
military forces are not near parity and have relatively 
"large" numbers of initial combatants, a deterministic model 
may be used to model the dynamics of combat without losing 
any essential information about the battle's outcome. From 
the research reported here, the author concludes that one 
may safely take "large" numbers of initial forces (in par- 
ticular, of tanks) , to mean units of battalion size or 
larger. Thus, it appears as though current U.S. Army deter- 
ministic models such as the Bonder/IUA model are to be 
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preferred for sensitivity analysis because of their compu- 
tational advantages, once they are "calibrated" from Monte 
Carlo simulation output (for example, from DYNTACS or 
CARMONETTE) . 
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APPENDIX A 



COMPARISON OF TWO METHODS OF OBTAINING THE PROBABILITY OF WINNING 

In this appendix it is shown that two apparently differ- 
ent methods of obtaining the win probabilities are really 
the same. In his doctoral thesis, Springall [Ref. 28] 
derived a recursive expression for the probability of 
winning in a homogeneous force, constant coefficient sto- 
chastic Lanchester model. His expressions and the develop- 
ment of them are quite different from those of Brown [Ref. 7]. 
Brown developed the familiar set of partial difference 
equations : 



P (m,n:r,s) 



A(r ,s) 



P (m, n : r-1, s) 



A(r, s) + B ( r , s ) 



B (r ,s) 



P(m,n:r,s-1) (66) 



A ( r , s ) + B (r , s) 



with boundary conditions 




otherwise 



n 



s 



(67) 



P(m ' n bp !r,n bp ) 




m = r 



otherwise 



( 68 ) 



P(m bp ,n:r,s) 



0 n > s; 



Ptm,!^ :r,s) 



0 m > r 



(69) 
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(m bp' n bp :r ' s) = 



= r 



■"bp 
otherwise 



and n. = s 
bp 



(70) 



where P (m,n bp :m o ,n o ) = probability X will win with m 
survivors. When this is summed from +1 to m Q , it yields 
the probability X will win: 



m 



X 



= Z 
m=: 



“b] 



+ 1 



P(m ' n bp :m o' n o ) 



(71) 



When the above equations are looked at carefully, what 
they are doing is explained by the following figure (Fig. 

74) . 

The probability that X will win with m survivors is 

simply the sum of the probabilities of all possible paths 

from ( m 0 / n 0 ) to (m,n bp ) with the probability of going from 

state (r , s) to state (r-l,s) equal to cffy ^ (± m m s ) and 

the probability of going from state (r,s) to state (r,s-l) 

equal to •%—. — . . And the probability X will win 

a t r , s ; + B(r,s) 

is simply the probability that X will win with m Q survivors, 
or m Q -l survivors, or ..., or +1 survivors. 

A similar analysis of Springall's equations will show 
that the two methods are identical. Springall concludes 
that 



P(m ' n bp :m o' n o ) = B(m,n bp +l)X(m,n bp +l) 



(72) 



n 






n o - 1 



r>o 



n bD+l 



n bp 



m bp 



+ 



• • • 



4 - . 



V 

+• 



+ 



m bp 



+r 



PlO 

ii 

Ai^.a^o). 



A(ni 0 ,no)+B(m 0 ,no) 

— J— «; — |- «s 1 



V 

+ 



+ 



V 

+ 



V 

+ 



-j- -c -f- 



m 0 -2~ 



“o -1 



Poi 

II 

-} ~ B(mo,no) 

A(m ,n 0 )+B(m 0 ,no) 



P02 



+• 



A(m 0 ,n 0 -1)+B(m 0 ,n 0 -1) 



+ 



n 



m. 



Method of Solving Recursive Relationship for the Probability 
of Winning for Brown's Equations 



FIGURE 74 
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where X(m, 11^+1) is derived from the following set of 
recursive relations; 



X(m o' n o } A(m ,n ) + B(m ,n ) 
0 0 o o 



B(m ,n+l)X(m ,n+l) 

X(m o' n) A (m , n) + B(m ,n) n bp < n < n o 



X(m,n Q ) 



A(m+l,n o ) X(m+l,n o ) 
A (m, n Q ) + B (m / n Q ) 



m < m < m 
np o 



X (m,n) 



_ A(m+l,n) X(m+l,n) + B (m,n+l) X (m, n+1) 
A (m,n) + B(m,n) 



Define 



P‘ (r,s) = (A ( r , s ) + B(r,s) )X(r,s) 



Then 



p, t m o' n o ) = 1 

B(m ,n+l) 

p,(m o' n) ~ A(m 0 ,n°l) + - B(m o ,nVl) p '‘ ro o- n+1 > 

A(m+l,n o ) 

P ' (m ' n o> ° Arm+r,np - + - B( - m Tr , i g - pl(l,>fl ' n o > 



( 73 ) 



( 74 ) 



( 75 ) 



( 76 ) 



( 77 ) 



( 78 ) 



( 79 ) 



( 80 ) 



145 



P ' (m, n) = 



A(m+1 , n) P ' (m+1 , n) 



B(m,n+l)P' (m,n+l) 



A(m+l,n) + B(m+l,n) A(m,n+1) + B(m,n+1) 



Therefore 



( 81 ) 



P (m. , n :m , n ) 

bp o o 



Atn^p+l,!!) P ' (m^p+1 , n) 
A(m bp +l,n) + BCit^p+l,!!) 



(82) 



p(m ' n bp :m o' n o ) 



A(m,n bp +l)P ' (m,n bp +l) 
A(m,n bp +1) + B(m,n bp +1) 



(83) 



What these equations say is, the probability of starting 

with (m ,n Q ) survivors is one. The probability of having 

A(m , n ) 

one more casualty from X is rr-i r — t-ttt \ and from 

B(m ,n ) o o o o 

Y is 1 ^ . . . The probability of having (m,n) 

A (m ,n ) + B (m ,n ) ^ ^ 

o o o o 

survivors is the probability of having (m+l,n) survivors 
times the probability of an additional casualty being from 

X, added to the probability of having (m,n+l) survivors 
times the probability of an additional casualty being from 

Y. This recursive relationship leads to the following 
figure (Figure 75) , which is the same as Figure 74. Thus, 
even though the sets of relationships appear to differ, 
they do yield the same results. This was further supported 
by a comparison of results from Lee and Wannasilpa [Ref. 20] 
and the results from the enclosed computer program, which 
was set up using the recursive equations of Springall. For 
the same parameters, the results were the same within 0.00001. 



n 




FIGURE 75 



The differences are assumed to have occurred because of 
roundoff errors within the computer system used. 

The program accuracy was further checked by 



m 



n 



o o 

E P(m,n :m ,n _) + E 
m=m b +1 bp o o 



n=n, +1 
bp 



P (m^pf n :m Q , n Q ) = 1 (84) 



All results were in the interval (0.9999692,1.0000000). 



147 



APPENDIX B 



NUMERICAL SOLUTION OF THE FORWARD KOLMOGOROV EQUATIONS 

Exact solutions for the time state probabilities for 
the square-law stochastic attrition process can be developed 
only for specific initial force levels. To the author's best 
knowledge, no solution for the time state probabilities for 
general initial force levels has been developed. Even if 
one were able to develop such a general solution, the result 
would be of no practical use due to its complexity. The 
solution for the linear-law stochastic attrition process 
was developed by Clark [Ref. 8], but it is too complex for 
practical use: So, solutions, if they could be developed, 

for more complex Lanches ter- type stochastic attrition pro- 
cesses, would most likely be of no practical use. Therefore, 
the most practical way to solve such equations is through 
numerical integration techniques. The method chosen for 
this thesis was the fourth order Runga-Kutta method [Ref. 21] . 
The computer program is set up to solve Lanchester-type 
attrition process models of the form 



dP (t ,m,n) 
dt 



A(m+l,n) P(t,m+l,n) + B (m, n+1) P ( t , m, n+1) 

- (A(m,n) + B (m,n) ) P (t,m,n) (85) 
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with 



P (0 ,m,n) = 



1 n = n and m = m 
o o 

0 otherwise 



( 86 ) 



P ( t , m, n) =0 m < in or n < 

' 9 o o 



(87) 



and boundary conditions 
dP ( t,m,rx ) 

dt “ = BI.,yl|P|t,.,n bp tl) 

dP(t,m ,n) 

St 22 Aln^p+l.nJPtt.^p+l.n) 



( 88 ) 



(89) 



However, with a few minor alterations, the program will 
solve the more general Lanches ter- type stochastic attrition 
processes with attrition rates of the form F(t,m,n) and 
G(t,m,n) . 

Numerical precision of the Runga-Kutta methods is 
dependent on the length of the time step chosen. For the 
time steps chosen, numerical precision appeared to be 
satisfactory since 

m n 

0 o 

1 £ P(t,m,n) was always in the interval 

m=m, n=n, 
op bp 

(0.9992,1.0003) . 
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Because the expected number of survivors was calculated 
from 



m(t) = E E mP(t,m,n) (90) 

m=in bp n=n bp 

it was felt the solutions were accurate enough. Additionally, 
because exact solutions exist for general initial force 
levels for P(t,m Q ,n) and P(t,m,n o ), they were used to 
check the accuracy of the approximate solutions . The 
approximate solutions never deviated from the exact solution 
by more than 0.00005. 
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APPENDIX C 



SCALING 

In a comparison of a stochastic model and the equiva- 
lent deterministic model, two classes of comparisons need 
to be considered. The first is the difference in results 
obtained from the two models. The second is to take these 
differences and see how they change with variations in 
parameters. One runs into difficulties when trying to make 
sense out of these comparisons when results are expressed 
in dimensional units (i.e.. What is the significance of 
the fact that at 150 minutes into a battle, the difference 
in survivors is 1.8 tanks?). 

One method of overcoming this problem is the technique 
of scaling. If comparisons are made when results are in 
dimensionless quantities, it is much easier to see if there 
are any significant differences. For example, one can 
easily compare two probability statements and note any 
difference between them. 

Scaling is a subset of a "science" called dimensional 
analysis, a technique often used in the physical sciences. 
For example, in the area of fluid dynamics, Reynolds deduced 
from some experiments that laminar flow broke down into 
turbulent flow at some critical velocity above that at 
which turbulent flow was restored to the laminar condition; 
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the former velocity being called the upper critical velocity, 
and the latter, the lower critical velocity. Reynolds was 
able to .generalize his conclusions by the introduction of 
a dimensionless term called the Reynold's number. He found 
that certain critical values of the Reynold's number defined 
the upper and lower critical velocities for all fluids 
flowing in all sizes of pipes, and thus deduced that single 
numbers define the limits of laminar and turbulent pipe 
flow for all fluids. Many other dimensionless numbers are 
used to generalize results — the Froude number, the Cauchy 
number, the Mach number, etc. Even in mathematics dimen- 
sionless quantities are used to generalize results — it and 
e, for example. (ir is the ratio between the circumference 
of a circle and its diameter, regardless of the units of 
measure employed.) 

Scaling involves the ratio of two quantities with the 
same dimensions. The dimension of a ratio of two quantities 
with the same dimensions is said to be dimensionless. The 
analysis in this paper is concerned with two classes of 
dimensions; time and quantity (number of casualties/survivors, 
and the differences in casualties/survivors) . Either model 
(the stochastic or deterministic) could have been chosen as 
the basis for forming ratios for the two classes of dimensions 
the deterministic model was used as a basis. Time was trans- 
formed into a dimensionless quantity by using the ratio t/t^, 
where tf is the time at which the deterministic battle would 
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end. It is easy to see that if 0 <_ t <_ t^, then the ratio 
is between zero and one. 

The transformation of the quantities requires some 
explanation of the rationale behind it. The transformation 
used was 

. _ , _ number of casualties at t 

ratio number of casualties at t^ for the deterministic mode 



It is easily seen that, for the deterministic model with 

0 <_ t <_ t^; 0 <_ A rat j_ 0 <_ 1 . A dimensionless quantity was 

desired that would express the "importance" of the number 

of casualties. A difference of one casualty out of forty 

is not as "important" as a difference of one casualty out 

of four. ^ ra tio ex P resses this "importance". For example, 

at t f , A . . = 0 for the deterministic model, while 

f ratio ' 

A ra tio 7*0' general, for the stochastic model. In fact, 



the difference in casualties at t. 



ratio 



number of casualties at t^ for the deterministic model' 



as the number of casualties at t^ for the deterministic 
model increases, A rat ^ Q decreases, giving the desired result. 
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APPENDIX D 



DERIVATION OF EXPRESSION FOR BIAS IN THE SQUARE LAW STOCHASTIC 
ATTRITION PROCESS WITH FIXED BREAKPOINT FORCE LEVELS 



From previous results it is known that the square-law 
stochastic attrition process is defined by: 

for m, +1 < m < m and n, +1 < n < n ; 

bp — — o bp — — o 



dP(t,m,n) 

dt 



anP(t,m+l,n) + bmP(t,m,n+l) 



(an+bm) P ( t ,m,n) 
(91) 



and on the boundaries of the state space. 



dP(t,m,n bp ) 

dt 



bmP (t ,m,n bp +l) 



(92) 



dP (t,m np/ n) 

at 




also, 



(93) 



P(t,m,n) - 0 



for m > m 

o 



and/or n > n 

o 



and 



P (0 ,m,n) 






1 




m = m and n 
o 

otherwise 



n 

o 
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Define 



A m P(t,m,n) = P ( t ,m+l ,n) - P(t,m,n) 



and 



A n P(t,m,n) = P(t,m,n+1) - P(t,m,n) 
so that (91) may be written as 



for 



m^p+1 < m < 



m and n, +1 < n < n 
o bp — — o 



dP ( t ,m,n) _ . . , 

dt _ ^ / rn f n) 



m 



bmA P(t,m,n) 
n ' ' 



First consider two jointly distributed discrete random 
variables, M and N with probability mass denoted as 
f^ N (m,n) . Further, assume that f^ N (m,n) = 0 for m < 
and n < n^ . Then the expected value of M, denoted as 
E(M) is given by 

00 oo 

E(M) = Z Z mf (m,n) 

M,N 

m=m, n=n, 
bp bp 



(94) 



^bp 



This may be expressed in terms of the conditional expectation 



00 



m=m, 

E (M :N=n) = =£ 



Z mf M,N (m,n) 



f N (n) 



as E(M) = e f N (n)E(M:N=n) 

n=n, 

bp 



Note that 



m 



o n 



dm 

dt 



[— E 
L dt L 
m : 



a 



Z mP ( t,m,n) ] 

nv n=n 
op bp 



(95) 



dm 

To compute which is given by equation (101) below, first 

sum equation (94) over all n, from n^ +1 to n Q to obtain 



n 



n 



dt n=n bD +1 

bp 



P ( t ,m, n) ] = a Z 



n=n,+l 

bp 



nA^P ( t,m,n) 



n 



o 



+ bm Z 



n=n, + 1 
bp 



A r P ( t,m,n) 



(96) 



Note that 



n 



n 



Z A P(t,m,n) = Z [P(t,m,n+1) - P(t,m,n)] 

n=n, +1 n n=n, +1 

bp bp 



= P(t,m,n o +l) - P(t,m,n b +1) 



or 
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Thus 



n 



n=n, +1 
bp 



A P(t,m,n) = 



n 



- P (t,m,n bp +l) 



n 



n 



dt 



[ E P ( t,m,n) ] = a [ E 



n=n b P +1 



n=n, +1 
bp 



nA^P (t / m /n) ] 



- bmP (t,m,n bp +l) 



(97) 



Adding (92) to (97) results in 



n 



■^[ E P (t,m,n) ] 

n=n. 

op 



n 



= a [ E nA P ( t,m,n) ] 



n=n, +1 
bp 



m 



(98) 



Multiply (98) by m and sum over m from m^tl to m Q , and 
add (9 3) multiplied by : 



m n 

, o o 

ft‘ E 

m=m. n=n, 
bp bp 



E mP ( t ,m,n) ] = 



dm 

dt 



m n 

o o 



n 



= a [ E 

m=irL +1 n=n, +1 



E mnA P(t,m,n)] + m, E anP(t,m, 

m, +1 n=n, +1 m ^ n=n ^ 

bp bp 



bp 



(99) 



Observe that 



+ l,n) 
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m 



n 



m=m bp +1 n=n bp +1 



Z mnA P(t,m,n) 

m 



m 



= Z 



n 



mA Z 
m 



nP (t ,m,n) 



m=m, +1 n=n, +1 

op bp 



m 



n 



n 



o o o 

= Z m[ Z nP(t,m+l,n) - Z nP(t,m,n)] 

m=rn bp +1 n=n bp +1 n=n ’ — + 1 



bp 



m n 

o o 



m 



n 



= Z 



m Z nP(t ,111+1,11) 

bp M 



111=111^^+1 n=n^_+l 



o o 

Z m Z nP(t,m,n) 

m=m, +1 n=n, +1 

Dp bp 



m +1 
o 



n 



Z (m'-l) Z 
P 

m 



m=m b _+2 



n=n, +l 
bp 



nP ( t , m ' , n ) 



n 



o o 

Z m Z nP(t,m,n) 

m=itL +1 n=n, +1 
Dp bp 



m +1 
o 



n 



Z (m'-l) Z nP(t,m',n) 

"V 1 n=n b P +1 

m n 



o o 

Z m Z nP(t,m,n) 

m=m b P +i n=n b P +i 



n 



m, Z nP(t,m +l,n) 

* p n=n, +1 np 

bp 
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nP (t ,m' ,n) 



m n 

o o 

= Z (m'-l) Z 

m=m^ +1 n=n L _+l 



:n bl 



m n 

o o 



Z m Z nP(t,m,n) 

m=m b P +i n=n b P +i 



n 



VJL ""t'V 1 '" 1 



n=n, +1 
bp 



m 



n 



= £ (m'-l-m) E nP(t,m,n) 

m=m bp +P n=n bp +1 

n 



- m z nP(t,m +l,n) 

np n. +1 np 

bp 



m a 

o o 



= Z 



nP (t ,m,n) 



m=m, +1 n=n, +1 
bp bp 



n 



m. Z nP (t ,nv +1 ,n) 

** n=n, +1 
bp 



Substituting (100) into (99) 



dm 

dt 



m n 

o o 



= -a Z E nP(t,m,n) 

ra=m bp +1 n=n b P +1 



n 



• a V n C +1 np(t '”bp +1 ' n) 

bp 
n 



+ am, £ nP(t,m +l,n) 

n=n, n P 

bp 



(100) 
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m 



n 



= -a 



m=m bp +:L n=n b P +i 



nP(t,m,n) + am^n^P <t /% p +l .%,> 



m 

o 

= -a E 

m=m bp +1 



n 

o 

E nP(t,m,n) 



n=n 



bp 



+ 



m 

o 



a E 
m=: 



"bi 



n P(t,in,n ) 
+ 1 bp bp 



+ am b P n bp p(t ' m bp +1 ' n bp> 



m 



n 



n 



= -a Z 

m=m b! 



E 

n=n. 



bp 



nP(t,m,n) + a E 
n=n 



nP(t,m, ,n) 



bp 



m 



+ a E n P ( t ,m,n ) 

+1 bp bp 



+ 



am b P n bp P<t ' m bp + 1 ' n bp ) 




-an + 



n 



a [ E nP(t,m. ,n) 



n=n. 



bp 



+ n bp ro bp P(t ' m bP + 1 ' n bp ) 



+ 




E P(t,m,n )] 

m=m b P +1 



( 101 ) 
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Equivalently , 



— m o 

|£ = -bS + b[J mP(t,m,n bp ) + r^n^P ( t ' m b p'"bp +1 > 
m ^bp 



+ 



“bp 



E P(t,m. ,n) ] 

n=n, +1 Dp 
bp 



(102) 
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non onnnonnooooonoooonoooooooooooooo 



COMPUTER PROGRAM 



SOLUTION TO STOCHASTIC LANCHESTER EQUATIONS 
USING FOURTH ORDER RUNGA-KUTTA APPROXIMATIONS 



LG=NUM6tR OF TIME INTERVALS NEEDED FOR COMPLETE 
SOLUTION TO DETERMINISTIC EQUIVALENT EQUATION 
MC=NUM6ER OF COMBATANTS X STARTS WITH 
NC=NUMBER OF COMBATANTS Y STARTS WITH 
H = SI Z E OF TIME INCREMENT 

A<X,Y)=* 'ATTRITION* • RATE FOR X IN TERMS OF X-l £ Y- 1 

BIX, Y)=* 'ATTRITION' • RATE FOR Y IN TERMS OF X-l 6 Y- 1 

FX = BREAK PO I NT FORCE RATIO FOR X 

FY=BREAKPOINT FORCE RATIO FOR Y 

MdP= I 1 — FXJ *MO 

NB F= ( 1-FYJ *NO 

DIMENSION KXINKXYJ, KYiNKXY) 

DIMENSION PROB (M0+1-M6P,N0+1-N8P) 

DIMENSION WK ( MO* 1- MBP » NO + 1-NBP , 3 ) 

DIMENSION DlLQ+1 ) ,DETERM(LU+1) 

DIMENSION Pi ( LO* 1 , MO* 11 , P2 ( LO* 1 , MO+1 J 
DIMENSION EST(LG*l), VAR(LQ*1) 

IF A PLOT OF THE TIME HISTORY OF THE X SIDE IS DESIRED 
SET I T = 0 . OTHERWISE SET I T=THE POINT IN TIME AT 
WHICH A THREE DIMENSIONAL PLOT OF THE JOINT 
PROBABILITY DISTRIBUTION OF X AND Y IS DESIRED 
(NOTE- TI ME= I T*H) . 



DIMENSION PI (313,41) , P2(313, 41) , EST(313) ,D(313) , 

1FI2J , S I ZE (i2 ) ,KX( 100), KYI 100) , 

2PFC6(41 r 4j.) ,WK(4i,41,3) , 

3VAR(b13) ,DETERM(313) 

INTEGER X , Y , U , Z 

REAL*8 TTL(12) /• CRAIG J.*,5*' • , • SMC 1462', 5*' 

REAL XTITLE/'TIME'/,YTITLE/' SUR V '/,ZTITLE/* '/ 

REAL Ki ,K2 ,K3 , K4 ,K1A, K2 A , K 3 A , K4A , K2B , K3 B , K4 B , K 1C , K2C , 
1K1B,K1D,K2 D,KbD,K4D,K1E,K2c,K3E,K4E,K1G,K2G,K3G,K4G, 

2 K4F , K3 C , K4C , KlF»K2F,K3r 

LOG! CAL*1 I DN( MO+l-MBP, NO* 1-NBP J 

LOGICAL* 1 IDN( 41 ,41 ) 

F ( 1)=0. 0 
F ( 2 ) =0 . 0 
SIZE( 1 ) =7. 0 
SI2E( 2) =9* 0 
LC=31 2 
MC=40 
N0=40 
NKXY-100 
MC 1=M0* 1 
LC1=L0* 1 
NC1=N0* 1 
I T= 1 0 
A 1 = 0.008 
61=0.004 
H=2.0 
FX = 1 .0 
F Y = 1 .0 

MB F= ( 1-FX I *M0 
NBP= ( 1-FY J *N0 
MS=MO+ 1— MB P 
NS=NO* 1-NBP 
M6P1=MSP*1 
NBP1 = NB P* 1 
M S1=MS- 1 

600 DO 13 L=1 , L01 
EST(L) =0.6 



162 



VAR(LI=0.0 
13 CONTINUE 

OC 400 N=1 » NS 
Y=NG+2-N 
DO jOQ M=1,MS 
X=M0+2-M 
DO 210 L =1 » L01 
IFO_.EO.iJ 00 TO 500 
I F ( X . EO . MO +1 ) GO TO 502 
I F ( X . E U • MB P+ 1 ) GO TO 504 
IF(Y. EQ.NO+1 JGO TO 506 
IFlY.EQ.NBPD) GO TO 507 

Kl=A<X«-l,Y.)*Pi(L-l,X*l)+B(X,Y+lJ*P2(L-ltX)-<A(X,Y) 

1 +8 (X , Y ) )*P1( L-ltX) 

K2=A(X+1, Y)*0.5* ( PlILtXO) + P1( L-1,X+1) ) ♦ B ( X , Y + 1 ) *0 .5* 
1(P2(L,X)+P2(L-1,X) )-(A(X,YD-6(X,YD*(Pl(L-l,X)+H*Kl 
2*0.5) 

K3=A(X*1, Y)*0.5*< P1(L,X<- 1) fr PI ( L- 1 , X* 1 ) ) * B ( X , Y+ 1 ) *0 . 5* 
ltP2(L,XD-P2(L-l,XD-(A(X,Y)+b(X,YD*lPi<L-l,X) + H*K2 
2*0. 5) 

K4=A(X+1,Y)*P1<L ,X+l)+B(X,Y*l)*P2(L,X)-< A<X,Y)*8(X,YD 
1*(P1(L-1»X) +H*K3 ) 

Pl(L»X) = Pl(L-i,X) + (H/6.0)*(K.l + 2. 0*K2 + 2. 0 *K3 + K4) 

GC TO 200 

500 IF(X.EQ.MO+1.AND.Y.EQ.NO+1)GO TO 501 
P1(L,X)=0.0 

GC TO 200 

501 F1(L,X)=1.0 
GO TO 200 

502 IFIY.NE .NO+1 )G0 TO 503 

P1<L,X) = EXP(-(A<X,YDB<X,YD*<L-1)*H) 

GC TO 200 

503 I F ( Y.EQ.NBPDJGO TO 509 

K 1 A=B ( X , Y*- 1 ) * P2 ( L— 1 » X ) - ( A ( X , Y ) *B ( X , Y ) ) * P 1 ( L-l , X ) 
K2A=e(X,YJ-l)*0.5*(P2(L,X)+P2(L-l ,X) ) 

1-(A(X, Y)+B(X,Y) )*( P1(L-1»X) *H*K 1 A*0 . 5 ) 
K3A=B(X.Y*1)*0.5*( P2(LtX)*P2(L-1 ,X) ) 
1-(A(X,yWB(X,Y))*( PI (L-l » X ) +H*K2A*0.5) 
K4A=B(X,YJ-1)*P2(L, X) 

1-(A(X,YDB(X,Y) )*( P1(L-1,X) +H*K3A) 

PI (L,X) =T*i ( L- 1 ,X)MH/6.0J*(KlA + 2. 0*K2A + 2 . 0 *K3A+K4A ) 
IF(Y. EQ. NO. ANO.L .EQ.2) GO TO 601 
GC TO 200 

601 E = ABS( P 1< 2,M01 )-( ( B ( MOD , NOD ) *N0/ A ( MO* 1 ,NODD* 

1 ( EXP ( ( A (MOD , NO D ) / NO ) *H)-i ,0)*P2( 2, M01) ) ) 

IFIE.LE.O. OODGO TO 200 
H=H*0. 5 
GO TO 600 

504 IF(Y.EQ.NBPD)GO TO 505 

K IF = A( XD,Y)*Pl(L-l,XD) 

K2F=A(XD,Y)*0.5*(P1(L,XD) +P1(L-1,XD) J 
K3F=A( XDtY)*0.5*(P1<L,X+1)+P1(L-1,XDD 
K4F = A(XD,.Y)*P1( L,XD) 

PI (L,X) =P1 4 L-l ,X ) + (H/6 ,0)*( K1F + 2 . 0 *K2F + 2 . 0 *K3F+K4F ) 

GO TO 200 

505 P 1 ( L i X ) =0. 0 
GC TO 200 

506 K1D=A( XD,Y)*Pl(L-l,XD)-(A(X,Y)+8(X,YD*Pi(L-l,X) 
K2D=A( XD,Y)*0.5*(Pi(L,X+l)+Pl(L-l,XDD- 

1 ( A (X t Y ) +B1X,Y) )*( Pl( L-l ,X) +H*K1D*0.5) 

K3D= A( XD* Y) *0.5*( Pi (L ,X+ 1) + Pi( L-l ,X+1D 
1-(A(X,Y)+B(X,YD*(PHL-1,X) +H*K2D*0.5) 

K4D=A( XD,Y)*P1(L,XD) 
l-( A(X, Y) +B( X, Y) ) *( Pi( L-l. X)+H*K3D) 

Pi (LiX) =Pi(L-l,X )DH/ 6-0)*(K.lD + 2.0*K2D + 2. 0 *K3D+K4D ) 
IF<X.EQ. MO. ANO.L .EC.2)G0 TO 602 
GO TO 200 

6 02 E = ABS ( P 1 ( 2-, MO) - ( (A(MOD,NOD)*MO/B(MOD,NO+1) )* 

1 <EXP< ( 6 (MOD, NOD) /MO) *H)-1 .0) *P1< 2, MOD )) 
IF(E.LE.0.001)GO TO 200 
H=H*0. 5 
GC TG 600 
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no 



507 K1G =B!X,Y+1)*P2! L-l ,X) 

K2G=B(X,YJ-1) *0.5*! P2!L,XH-P2!L-1 ,X) ) 

K3G=B( X,Ym*0.5*!P2(L,X)+P2!L-l,X)) 

K4G=6! X,Y.«-1)*P2! L,X) 

P1(L»X) = PJL IL-i ,X)«- !H/6.0)*!KlG + 2.0*K2G + 2 .0*K3G*K4G) 
GC TO 200 

509 K1G =6<X,Y*1)*P2(L-1,X) 

K2G=B( X,Y*l)*0.5*l P2 ! L , X )+P 2 ! L- 1 ,X) ) 
K3G=B<X,Y*1)*0.5*! P2!L,X)+P2!L-1,XJ ) 

K4G=B! A, Yj-1)*P2! L,X) 

PI <L,X) =Pl!L-i ,X)MH/6.0)*< K1G+2.0*K2G + 2.0*K3G<-K4G) 
200 CONTINUE 
210 CONTINUE 
300 CONTINUE 

DO 1 I = M8Plt M01 

IF! IT.EC.OIGO TO 50 

PRCB( I-MBP.Y-NBP )=P1UT , I ) 

50 00 1 J = 1 » LOl 
P 2! J » I ) =P1 1 J » I ) 

VAR! J)=VAR! J|t( I-1.0)*!I-1.0)*P1! J,I) 

EST( J)=ESJ!J)MI-1.0)*PUJ, I) 
l CONTINUE 
Z=Y~1 

I F ( I T . GT .0 )G0 TO 901 
IF(N01-MBP.GT.15)GO TO 901 

WRITE!6,100) Z, ( I , I=MBP, MO) 



901 DC 99 J = 1 * LOl 
D( J)=( J-1J-H 

IF ! IT..GT .0 )GG TO 99 

IFIMOl-MBP.GT .15 ) GO TO 99 

WRITE (6,101)0! J) f< PlUiIl ,I=MBP1»M01) 

99 CONTINUE 
400 CONTINUE 

TI = ( IT-1).*H 
WPITE(6,1G5)H 
WRITE! 6,103) 

A2=SQRT ( Al*8 1 ) 

A3-SQRT ( Ai /31 ) 

A4= 1 . 0/ A3 
A5=MG- NG*A3 
A6=MG+ NO* A3 
DO 2 J= 1 , LG1 
LC = J 

F1=EXP!A2*D( J) ) * A5 
F2=EXP!-A2*D( J) ) * A6 
DETERM(J)?0.5*(F1+F2) 

VAR! J) = VA-R! J)-E$T(J)*EST(J) 
WRITE(6,102)D(J) , E ST ( J ) , DE TERM ( J ) , VAR! J) 
IF(DETERM! J) .LE.M3P)G0 TO 850 
YFORCE=0.5*A4*(F2-F1) 

IF!YFORCE-LE.NBP)GO TO 851 

2 CONTINUE 
GO TO 700 

850 Iw=l 

GO TO 700 

851 IW=0 

GO TO 700 

700 IF! I T.EQ.OJGO TO 900 
WRITE! 6 , 1 Oc ) T I 
WRITE! 6,107) ( NT , NT=NB P , NO ) 

MJ=«Ol-MdP 
NJ=N01-N8P 
DC 110 1=1, MJ 
MT-I+MBP-1 

110 WRITE(6,104)MT,( PRCB!I ,J) , J=1,NJ) 

900 IF ! IT-GT.OIGU TO 902 
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DMIN=DETERM(LCI 
TF1=0.25*D*LC)/H 
TF2=2.0*TF 1 
TF3=3.0*TF1 

WRITE*6,108)M0,N0,MBP,N6P,D(LC),TF1,TF2,TF3 
I F ( I W . EU .1 ) GO TO 852 
WRITE* 6, 109) 

GO' TO 8 53 

852 WRI T E ( 6 1 11 1 ) 

853 IFIEST(LC) .LT. DETERM* LC) ) GO TO 270 
DM 1= MO— 0M1 N 

DO 903 J = 1 » LC 

DETERM* J) = l * DETERM* JJ-DMIN) /DM1 ) *100.0 
D(J)=D(Jl#100.Q/D*LC) 

903 CONTINUE 
CALL PLOTS 

CALL SCALE* DETERM, LC,1 , 9 .0 , 1 . 0 , YM IN , DY ) 

CALL SCAL5*D,LC,i , 6. 0,1.0, XMIN,DXJ 

CALL AX IS* 0. 0,0.0, XT I TLE, -4,6. 0 , 0 .0, XMI N , DX ) 

CALL AXIS* 0.0,0.0,YTITLE,4,9.0,90.0,YMIN,DY) 

CALL LINE* D, DETERM, LC, 1,1) 

DO 14 1=1, LC 

EST * I ) = l ( * EST* I)-DMIN)/DM1)#100.0-YMIN) /DY 

14 CONTINUE 

CALL LINE*D,EST,LC,1,1) 

CALL PLOTE 
GO TO 271 

270 DM1 = MO- DM! N 
DO 950 J=1 , LC 

DETER Ml J i = < * DETERM* D-DMIN) /DM1 1*100.0 
D*J)=D( J)£100.0/D(LC) 

950 CONTINUE 
CALL PLOTS 

CALL SCALE (DETERM, LC,1 ,9.0 . 1 .0, YMIN.DY ) 

CALL SCALE*D,LC,1»6.0,1.0,XMIN»DX) 

CALL AXIS*0.0,0.0,2TITLE,-4,6.0,0.0,XMIN,DX) 

CALL AXIS* 0.0,0. 0, YTITLE,4 ,9 .0 , 90 .0 , 0 . 0 , DY ) 

CALL LI Nc*D, DETERM, LC, 1,1) 

DO 30 1 = 1 » LC 

EST(I)=**(EST( I ) — DM I N) /DM1 ) * 100 . 0— YMI N) /DY 
30 CONTINUE 

CALL LINE*D,EST,LC,1,1) 

CALL PLOTE 

271 IF* IT.EG.OJGO TO 999 
902 00 15 K = 1 ,.MS 

15 EST*K)=K-1 
DO 16 K= 1 , NS 

16 OETERM* K)fK- 1 
AA=15 . 0 
BB=30 . 0 
PMAX=PROB.< 1 , l) 

DC 18 1=1, MS 
DO 18 J = 1 , NS 

IF(PR03(It J) .GT .PMAX) P MAX = PR OB ( I ,J) 

18 CONTINUE 
CCNST=8.0yPMAX 
DO 19 1=1, MS 
DO 19 J = 1 , NS 

PROB ( I , J)=CGNST*PR06*I , J) 

19 CONTINUE 

CALL PLT3D1* EST, MS, DETERM, NS, PROB, AA,BB, F,TTL, SI ZE, 
1WK* IDN ,KX,.KY ,NKXY, 0) 

999 CONTINUE 

100 FORMAT* IX////, 3X, ' THE FOLLOWING ARE FOR • , I 3, • SURVI', 
l'VORS ON THE Y 6 I D E ' // / , 20 X , • NO . OF X SURVIVORS'//, 
26X,'TI ME* , 1517////) 

101 FORMAT (6X,F7.3»15F7.4) 

102 FORMAT* l6X.,F7.3,15X,F11.7,15X,F11.7,18X,F11.7) 

103 FORMAT * 1 5X , * TIME ' , 1 5X , • E XPE CT ED NO. OF SURVIVORS', 

1 5X, ' DETERMINISTIC EGU I VALENT ', 5X , ' VARI ANCE* /) 

104 FORMAT ( 13X , 1 7 , 1 5 F7 .4// ) 
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105 FORMAT! IX///, 5X, 'H IS ',Fll.7////J 

106 FORMAT! IX-,/////, 3X , 'MATRIX OF PROBABILITIES OF NUMBER' 
1' GF SURVIVORS AT TIME ',F11.7//i 

107 FORMAT! 19X,15I7//J 

108 FORMAT! IX///, 5X, • INITIAL FORCE SUE FOR X IS ',13/, 
15X, 'INITIAL FORCE SIZE FOR Y IS • , 13/ , 5X , • 6REAKP0 I NT • , 
2' FORCE LEVEL FOR X IS • ,I3/»5X, 'BREAKPOINT FORCE 

3 ' LEV EL FOR Y IS • , 1 3 / , 5X , • BA TT LE WILL TERMINATE AT ', 

4 ' TIME • ,FU.7///,5X,'TO COMPUTE 3-D GRAPHS AT 25, ', 

5 ' 50 , AND 75* TIMES INTO THE BATTLE, IT SHOULD BE'/.lOX, 
6 ' FOR 2 5%' ,.F 1 1 . 7 / , 1 OX » ' F OR 50% ' , F 1 1 . 7/ , 10X , • FOR 7 5% * , FI 

109 FORMAT ( 5X,. 'THE WINNER IS X'////J 
111 FORMAT (5X-, 'THE WINNER IS Y'////) 

END 

FUNCTION A! I ,J) 

A=0.00d*( J-l .0) 

RETURN 

END 

FUNCTION B ( I » J ) 

B=0. 004*1 1-1 .0) 

RETURN 

END 



ooooooo 



COMPUTER PROGRAM 



SOLUTICN TO SPRI NGALLS RECURSIVE SOLUTION FOR THE 
PROBABILITY ONE SIOE WILL WIN 



DIMENSION XI I50 t 40) ,NS(5) ,RS( 6) 



NS ( 21 =5 
NS( 3 ) = 1 0 
NSI4) =20 
NS (5) =40 
RSI 11=0.23 
RS(2)=0.5 
RSI ^J=l .0 
RSI4)=2.0 
RS ( 5) =3 . 0 
RS(6i=4.0 
RSI 7)=5.0 
RS I 8) =6 . 0 
FX=1 .2 
FY=1 .2 
DO B Ki = l-» 5 
FX=FX-0.2 
F Y=FY-0 . 2 
DO 6 L= I » 5 
NC=NS ( L ) 

DO 7 K= 1 1 8 
R=R S I K ) 

SM=2 . 0*NO*NO* I 1 .0-1 1.0-FYJ*! 1 . O-FY ) )*R/ I 1 . 0- 1 1 . 0- FX ) 
1*1 1.0-FX) ) 

MC=SQRT( SMI 
MBP= ( 1 . Q-FX) *M0 
- NB P= I l. 0-FY)*N0 
MSPl=MBP+i 
N6P1=NB P + 1 
MDIFF=MO-MBP 
NDI FF=NG-NBP 
MCOUNT=MDI FF-1 
NCCUNT=NDI FF-1 

X(MD1FF»NDIFF)=1 .0/(AIMO»NO)+BIMO»NOJ ) 

DO 1 J=1 i NCOUNT 
NREAL=NO-J 
NSTOK t = ND I FF-J 

X I MD IFF , NSTORE ) = I 6 1 MO t NRE AL + 1 ) *X IMDI FF » NSTORE + 1 J )/ 

1 I A I MO, N REAL) +B I M0» NRtAL ) ) 

1 CONTINUE 

DO 2 I = 1 » MCOUNT 
MREAL=MO-I 
MSTORE=MDI FF-I 

X(MSTORE.NOIFF)=( AIMREAL + 1 , NO ) *X I MSTORE + 1 , NO IFF ) J / 

1 I A IMRE AL, NO +0 1 M RE AL , NO ) ) 

DO 3 J=l, NCOUNT 
NREAL=NO-U 

NS TORE= NDI FF-J 

X I MS TORE ,N STORE) =( AIMREAL+ 1 , NRE AL) *X I MSTORE+ 1 , NSTORE ) 
1 +B(MREAL,NREAL + 1 ) *X { MS TORE » NSTORE+ 1) )/ 

2 I A IMRE AL » NREAL ) + B I MRE AL , NRE AL ) ) 

3 CONTINUE 

2 CONTINUE 
PMWI N=0 . 0 
PNW I N=0 . 0 

DO 4 I=1,MDIFF 

PMWIN = PMWIN + B( I +MBP , NBP 1 ) *X I I , 1) 

4 CONTINUE 

DO 5 J = 1 » ND IFF 

PNWIN=PNWIN+A(MBP1 , J+NBP)*XI 1, J) 

5 CONTINUE 
SUM=PMW I N+PNWI N 

S=MO*MO*( i.-( l.-FX)*( l.-FX) ) /I 2 .*( 1 . - { 1 . -FY ) * I 1 . -F Y) ) 
1*N0*NU) 
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WR I T E ( 6 i 100J MO,MBP ,N0,N8P ,FX,FY 
WRITE(6, 101) PMWIN, PNWIN, SUM f S 

7 CONTINUE 
6 CONTINUE 

8 CONTINUE 

100 FORMAT ( 5X, *MO=» , 1 4 , 5X , * MB P= * , 14, 5X, *N0= ' , 14, 5X, • N3P 
1 1 4 » 5 X , * FX=* ,F6.2*5X» , FY=* » F6 .2/ I 

101 FORMAT (5X, *PMHIN=* ,F11.7,5X, • PNWIN=* ,Fll.7»5X» *SUM= 
lFil.7 f 5X»* R=» »F1 1. 7/// J 

ENO 

FUNCTION All t J ) 

A=C.G08*J 

RETURN 

ENO 

FUNCTION B(I,J) 

B=0. 004-1 

RETURN 

ENO 



BIBLIOGRAPHY 



1. Adkins, Richard C. f Analysis of Unit Breakpoints in Land 
Combat, M.S. Thesis, Naval Postgraduate School, March 
1975. 

2. Bach, R.E., Dolansky, L., and Stubbs, H.L., "Some Recent 
Contributions to the Lanchester Theory of Combat," 
Operations Research , v. 10, p. 314-326, 1962. 

3. Bishop, Albert B. , and Clark, Gordon M. , The Tank Weapons 
System , Systems Research Group, Department of Industrial 
Engineering, The Ohio State University, 30 June 1969. 

4. Bonder, Seth, "An Overview of Land Battle Modeling in 
the U . S . " , Proceedings, Thirteenth Annual U.S. Army 
Operations Research Symposium , p. 73-88, 29 October - 
1 November 1974. 

5. Bostick, Steven P. , Brandi, Francis X., Burnham, C. Alan, 
and Hurt, James J. , "The Interface Between DYNTACS-X 

and Bonder I.U.A.", Proceedings, Thirteenth Annual U.S. 
Army Operations Research Symposium , p"! 494-502 , 29 
October - 1 November 1974. 

6. Brooks, F. , "The Stochastic Properties of Large Scale 
Battle Models", Operations Research, v. 13, p. 1-17, 

1965. 

7. Brown, Richard H., "Theory of Combat: The Probability 

of Winning", Operations Research, v. 11, p. 418-425, 

1963. 

8. Clark, Gordon M. , The Combat Analysis Model , PH.D. 
Dissertation, The Ohio State University, 1969. 

9. Dolansky, L. , "Present State of the Lanchester Theory 

of Combat", Operations Research , v. 12, p. 344-358, 1964. 

10. Engineer Strategic Studies Group, Bibliography of 
Publications , 1972. 

11. General Research Corporation, Operations Analysis 
Division, Gaming and Simulation Department, A Hierarchy 
of Combat Analysis Models . 

12. Hannah, William P., Further Comparisons of Stochastic 
and Deterministic Models for the Optimal Control of~ 
Lanchester-Type Attrition Processes , M.S. Thesis, Naval 
Postgraduate School, March 1975. 



169 






\ 



\ 




13. 



Helmbold, R. , Decision in Battle: Breakpoint Hypothesis 

and Engagement Termination Data ^ The RAND Corporation 
Report R-772-PR, June 1971. 

14. Isbell, J.R. , and Marlow, W.H., "Attrition Games", 

Naval Research Logistics Quarterly, v. 3, p. 71-94, 

1956. 

15. Jain, J.C., and Nagabhushanam, A., "A Two-State Markovian 
Correlated Combat", Operations Research, v. 22, 

p. 440-444, 1974. 

16. Johnson, Patricia C., Plotting Package for NPS IBM 360/67 , 
Naval Postgraduate School TN No. 02li-03, January 1974. 

17. Kisi, T., and Hirose, T. , "Winning Probability in an 
Ambush Engagement", Operations Research, v. 14, 

p. 1137-1138, 1966. 

18. Koopman, Bernard O. , "A Study of the Logical Basis of 
Combat Simulation", Operations Research, v. 18, 

p. 855-882, 1970. 

IS. Lanchester, F.W. , "Aircraft in Warfare: The Dawn of the 

Fourth Arm-No. V. , The Principle of Concentration", 
Engineering , v. 98, p. 422-423, 1914 (reprinted in 
World of Mathematics , y. IV, p. 2138-2148, J. Newman, 
ed . , Simon and Schuster, New York, 1956). 

20. Lee, Woo Young, and Wannasilpa, Amnuay , Comparison of a 
Deterministic and a Stochastic Model for the Probability 
of Winning in a Two-Sided Combat Situation , M. S . Thesis , 
Naval Postgraduate School, September 1972. 

21. McCracken, Daniel D. , and Dorn, William S., Numerical 
Methods and FORTRAN Programming , John Wiley and Sons, 

Inc. , 1964 . 

22 . Morse, P.M., and Kimball, G.E., Methods of Operations 
Research , Chapman and Hall, London, 1950. 

23. Naddor, Eliezer, "Dimensions in Operations Research", 
Operations Research , v. 14, p. 508-514, 1966. 

24. Powers, Robert L. , and Taylor, James G. , Comparison of a 
Deterministic and a Stochastic Formulation for the Optimal 
Control of a Lanchester-Type Attrition Process , Paper 
presented at the International Symposium on Applications 
of Computers and Operations Research to Problems of 
World Concern, Washington , D. C. , 20-21 August 1973. 



170 



25 . 



Raney, Sharon D. , PLT3D1: Three Dimensional Isometric or 

Perspective Off-Line Plotting Subprogram with Hidden Line 
Elimination , Naval Postgraduate School Technical Memorandum, 
February 1974. 

26. Smith, David G. , "The Probability Distribution of the 
Number of Survivors in a Two-Sided Combat Situation", 
Operational Research Quarterly, v. 16, No. 4, p. 429-437, 
1965. 

27. Snow, R.N. , Contributions to Lanchester Attrition Theory , 

The RAND Corporation Report RA-15078, April 1948. 

28. Springall, Anthony, Contributions to Lanchester Combat 
Theory , PH. D. Thesis - ) Virginia Polytechnic Institute, 

March 1968. 

29. Taylor, James G. , Lanchester-Type Models of Warfare , 

Tutorial delivered at the Thirty-fifth Military Operations 
Research Symposium, U.S. Naval Academy, Annapolis, 

Maryland, 1-3 July 1975. 

30. Taylor, James G. , "Solving Lanchester-Type Equations for 
'Modern Warfare' with Variable Coefficients", Operations 
Research, v. 22, p. 756-770, 1974. 

31. Taylor, James G. , "Survey on the Optimal Control of 
Lanchester-Type Attrition Processes", presented at The 
Symposium on the State-of-the-Art of Mathematics in 
Combat Models (also Technical Report NPS55Tw70431, 

Naval Postgraduate School, Monterey, California, March 
1974) . 

32. Taylor, James G., and Parry, Samuel H., "Force Ratio 
Considerations for Some Lanchester-Type Models of 
Warfare" , Operations Research , v. 23, p. 522-533, 1975. 

33. U.S. Army Concepts Analysis Agency, Methodology and 
Resource Directorate, Tabulation of Models of Interest 
to USACAA , June 1975. 

34. Vennard, John K., Elementary Fluid Mechanics , 3rd ed. , 

John Wiley and SoniT) Inc. , 1954. 

35. Weiss, George H. , "A Comparison of a Deterministic and 
a Stochastic Model for Interaction between Antagonistic 
Species", Biometrics , v. 17, p. 595-602, 1963. 

36. Weiss, Herbert K., "Lanchester-Type Models of Warfare", 
Proceedings, First International Conference on Operational 
Research , p. 82-99, Oxford, September 1957. 



171 



37. Willard, D. , Lanchester as Force in History: An Analysis 

of Land Battles of the Years 1618-1905 , Research Analysis 
Corporation Report RAC-TP-75, November 1962. 



172 



INITIAL DISTRIBUTION LIST 



No. 



1. Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 

2. Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 

3. Department Chairman, Code 55 
Department of Operations Research and 

Administrative Sciences 
Naval Postgraduate School 
Monterey, California 93940 

4. Associate Professor James G. Taylor, Code 55Tw 
Department of Operations Research and 

Administrative Sciences 
Naval Postgraduate School 
Monterey, California 93940 

5. Cpt. James D. Craig 
1321 West College 
Independence, Missouri 64050 

6. U.S. Army Military Personnel Center 

ATTN: DAPC-OPD-PD-CS 

200 Stovall Street 
Alexandria, Virginia 22232 

7. Chief of Naval Personnel 

ATTN: Pers lib 

Department of the Navy 
Washington, D.C. 20370 

8. Associate Professor Alan R. Washburn 
Department of Operations Research and 

Administrative Sciences 
Naval Postgraduate School 
Monterey, California 93940 

9. Dr. Wilbur B. Payne 

Deputy Undersecretary of the Army 

(Operations Research) 

2E621 Pentagon 
Washington, D.C. 20310 



Copies 

2 

2 

2 

1 

1 

1 

1 

1 

1 



173 



1 



10. Associate Professor Gordon M. Clark 
Department of Industrial Engineering 
The Ohio State University 
Columbus, Ohio 43210 



174 




The effect of uncer- 
tainty on Lanches*c r - 
type equations of com- 
bat. 



Thes i s 

C782< 

c.l 



thesC7824 

The effect of uncertainty on Lanchester 




3 2768 001 02441 7 
DUDLEY KNOX LIBRARY 





